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ABSTRACT 

This report is a user's manual for the F[ail-Safe Abort System TEMP erature 
Analysis Program, FASTEMP., This program was used to analyze fail-safe abort sys- 
tems for an actively cooled hypersonic aircraft (Contract NASA-Langley , NASl- 13631) . 
FASTEMP analyzes the steady state or transient temperature response of a thermal 
model defined in rectangular, cylindrical, conical and/or spherical coordinate 
system. FASTEMP provides the user with a large selection of subroutines for heat 
transfer calculations. The various modes of heat transfer available from these 
subroutines are: 

o Heat storage 
Q Conduction 
o Radiation 

o Heat addition or generation 
o Convection 
o Fluid flow 

These modes of heat transfer can be simulated by the program in the solution 
of one, two, or three dimensional heat transfer problems. The modes may be modeled 
by rectaxigular , cylindrical, conical and/or spherical geometries. Any combination 
of modes and geometries is allowed for analysis of either transient response or 
steady state temperatures. 

The program obtains its solution using the backward-difference method. The 
maximum number of nodes that may be solved depends upon the thermal model and the 
solution method used with an absolute maximum of 9999 nodes. 

The program is written in a combination of Fortran IV and Assembler language, 
COMPASS, for the CDC Series 6000 and CYBER Series Computers. Minimum core require- 
ments are 47000 core locations on the Scope operating system, 
o 

This report contains general information on the program and is designed to aid 
the engineer in setting up a problem and obtaining solutions with a fast turnaround. 
Program structure, solution techniques, input data and program output are discussed. 
A sample problem is presented and program usage is discussed. The program input/ 
output is designed for utilization of customary engineering units. 
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Symbol 

A 

CP 

C 

P 

F 

h 

k 

n 

0 

Q 

AT 

T 

T 

m 

V 

AX 

e 

P 

0 

b 


NOMENCLATURE 

Definition 

Area, in^ 

CDC Central Processor 

Specific heat at Constant Pressure, Btu/lb-°R 
Radiation View Factor 

Heat Transfer Coefficient, Btu/hr-f t^-°R 

Thermal Conductivity, Btu/hr-ft-°R 

Mass Flow Rate, Ib/hr 

Stands for alphabetic letter 0 

Heat Flux, Btu/hr 

Temperature Difference, ®R 

Temperature, °R 

Boundary Temperature, °R 

Volume, in^ 

Conduction Path Length, in 
Emissivity 

3 

Density, Ib/ft 

Stef an-Boltzmann constant, 0.1714 x 10 , Btu/hr-ft -°R^ 

Superscript, means leave a blank in given card column 
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1. INTRODUCTION 

This user's manual for the FASTEMP Heat Transfer Computer Program contains 
general information on program usage. The heat transfer program consists of a 
Thermal Analyzer (control subroutine) and user called subroutines, which are des- 
cribed in this user's manual. The computer code was written in Fortran IV and 
Assembler language (COMPASS) for the CDC 6000 and CYBER series computers. The 
user called subroutines are described in Appendix A. ■ 

This program was developed to analyze transient and steady state 
heat transfer problems. A problem is set up by defining a model which is made 
up of a finite number of lumped elemental volumes or nodes. Nodes are connected 
with the appropriate heat transfer terms. The user may model a very complex 
problem and is only limited by the amount of core available and solution accuracy. 
(Solution accuracy is discussed in Section 7.) The user may set up one, two or 
three dimensional problems using the following modes of heat transfer: 
o heat storage 
o conduction 
o radiation 
o convection 
o heat flux 
o fluid flow 

The modes of heat transfer may be modeled with rectangular, cylindrical, conical 
and/or spherical geometries. 

The program obtains a solution for temperatures using the backward finite 
difference method for the heat balance equations. This method results in a set of 
simultaneous equations which are solved. Several matrix solutions are available 
which optimize core or computer time. The maximum number of nodes (equations) in 
the thermal model depends on the selection of heat transfer modes and the matrix 
solution chosen as well as available computer core. The absolute maximum number 
of nodes which may be run is 9999 nodes. This limit is imposed by programming 
constraints using a four column field to specify node numbers. 

The primary output of the computer program is a tabulation of the tempera- 
ture of each node at problem times specified by the user. Additional output is 
printed for some heat transfer subroutines as described in Appendix A. 

This user's manual is based on the program write-ups presented in References 
(1) and (2). The program reported herein and the referenced program are data 
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compatible although operational on different computers for over nine years. Refer- 
ence (2) program contains additional subroutines which allow building thermal 
models using phase change, aeroheating or orbital flux calculations, and other 
analysis oriented subroutines. 
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2. PROGRAM STRUCTURE 

FASTEMP contains a thermal analyzer control subroutine which interfaces between 
the user and the computer system. The user initiates the job by compiling a sub- 
routine (EQDAT) which calls heat transfer mode subroutines. The compiled EQDAT 
subroutine is loaded with the thermal analyzer controller and the library of user 
called subroutines and execution is initiated. The user then supplies data for the 
called subroutines and a model is analyzed by the program. A heat transfer problem 
may be analyzed using a full set of data, called a basic case, or a partial set of 
data which are changes to a previously run case. (Input for basic and change cases 
is discussed in Section 4.) 

Job control language (JCL) cards are also required for each job to load and 
execute the program. The JCL is a function of the computer system and is not covered 
in this report. 

After loading, the execution of the main program is initiated. The control 
subroutine calls the subroutine EQDAT which is required and compiled for each job. 
This subroutine in turn calls the subroutines needed to describe the model to be 
analyzed. The call establishes the order in which the subroutine data must be input 
to the program. An example of Subroutine EQDAT is shown below: 

SUBROUTINE EQDAT 
CALL SRA 
CALL DRA 
« 

. Subroutines called by user 

CALL RRA ' 

RETURN 

END 

In developing an effective tool useful for analysis of both small models (less 
than 20 nodes) and large models (over 100 nodes), the program was structured into 
input, setup and computing sections. 

These are designed as the M (minus), the ^ (zero), and the _F (plus) phases 
of the program as denoted by the sign of the control key. The M phase preprocesses 
the input data cards for a case. Error messages are written if any data cards are 
missing (if required), unrecognizable or in the wrong order. If tlie case being 
processed is a change case, the cards will be merged with those of the case being 
changed. The program has two levels of change case capability. Any number of 
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change cases may be run which change the basic case. A change case may also be 
run which changes the most recent change case. The ^ phase reads the numeric 
input values, performs error checks, computes constants using geometrical data and 
stores the data for use by the ^ phase. The ^ phase computes the heat flux terms 
for each node, solves the heat balance matrix for node temperatures, and controls 
the output of the program. 

Several computer clocking points are provided in the program. The computer 
-time at the beginning of the M phase will be printed with the message time 
expended = XX CP (SEC). The computer time at the end of the Z phase will be 
printed- with the same message. At the end of a case, the clocking printed is 
the amount of time for the ^ phase (time step calculations). All clocking times 
are computer central processor times and are printed in seconds. 

The Thermal Analyzer control subroutine uses a dynamic core allocation sub- 
routine which allows the core requirement for the program to vary in size as a 
function of the size of the model being analyzed. The program prints the number 
of locations used by the subroutines and that remaining of the core size specified 
for the job at several stages in the running of a case. The printed messages are 
self explanatory. 

The ^ phase is used to compute temperatures from data supplied by the user. 
The analysis is a time transient calculation at intervals composed of a user sup- 
plied time step schedule. A time step is processed in the following manner. The 
£ phase calls Subroutine TIMST to calculate the time step interval from user data. 
Then Subroutine EQDAT is called. Subroutine EQDAT calls the heat transfer method 
subroutines which compute and store the matrix coefficients for each heat transfer 
term. Properties which are a function of node temperature are evaluated at the 
known temperatures at the beginning of a time step. Values which are a function 
of problem time are evaluated for backward finite differences at the time at the 
end of the time step. Subroutine SOLVE is called after the matrix coefficients 
have been computed and stored by Subroutine EQDAT. Subroutine SOLVE calls the 
matrix solution subroutine (MSOLV), and the print subroutine (PRINT). Subroutine 
MSOLV is called to compute the node temperatures and Subroutine PRINT is called 
to output the results (if requested). The computed temperatures are stored as the 
known temperatures for use during the next time step, the matrix locations are set 
to zero, the problem time is updated, and the above computations are repeated. 
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3 . SOLUTION METHODS 


A heat transfer problem is analyzed by defining a model consisting 
of a finite number of nodes. The nodes or elemental volumes are considered 
to be isothermal at any given time during the analysis. A heat balance 
is written for each node using the heat transfer subroutines which are explained 
in Appendix A. Most of the terms in the heat balance are heat flux terms, 
however, some are energy terms, e.g. fluid transport terms (m AT). In this 
report the term heat balance includes heat flux terms and energy terms. The heat 
balance for a node is not restricted to any number of terms for a given mode of 
heat transfer, and the heat balance for a node may be written as: 

ZQ stored = ZQ conduction + ZQ convection + ZQ radiation + ZQ misc. 

All terms in the equations have the units of Btu/hr. Program temperatures are 
in degrees Rankine internally although input of initial temperatures and output 
of computed temperatures may be in either Rankine or Fahrenheit at the users 
option. A consistent set of units is used for dimensions, e.g., lengths are 
input in inches and areas in square inches except where otherwise specified. 

The thermal analyzer uses heat balance equations formulated for each node 
using the backward finite difference method (Reference (3)), A set of simultaneous 
equations results with coefficients stored in matrix form. The user may optionally 
select one of several different matrix solutions (called METHODS) to optimize the 
core requirements or solution time. 

3.1 BACKWARD DIFFERENCE SOLUTION (METHODS 3 THROUGH 8) 

The backward difference heat balance equation expresses all heat fluxes 
in terms of the temperatures and boundary conditions at the end of a time step. 

The general form of the backward difference heat balance for one node is: 


V (T' - T°) 


= M (X' - T') + hA (T' - T') + OeFA (T'^ 
AX m m m 


where T' = unknown node temperature at end of time step 

T° = known node temperature at beginning of time step 

T' = known boundary or unknown surrounding node temperatures at 
m 

the end of time step 


This expression is used to represent a set of simultaneous equations, linear in 

4 

in T with the exception of the T radiative heat flux terms. An approximate 
linearized term is substituted for radiative terms to allow the simultaneous solu- 
tion of equations which are linear in T. (The approximation used for linearization 
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is presented in Appendix B.) The backward difference equation solution is always 
stable. However, increasing the time step increases truncation errors in the 
solution and the user must tailor his time step to the boundary conditions. A 
general discussion of this aspect is given in Reference (3) . 

Thermal models produce matrices which can vary widely with respect to the 
connection of equations represented by the matrix coefficients. The thermal 
analyzer has several matrix solution methods to allow the user to specify the most 
efficient solution technique for a given problem. The efficiency of the specialized 
solutions results from reductions in computer core storage and/or running time. A 
brief description of methods 3, A, 6 and 8 supplied in this program is given below. 

Method 3 - This method uses a full square matrix. Method 3 may be 
used for any backward difference problem, however, it 
will waste computer core space for many problems. 

Method 4 - This method uses a tridiagonal matrix which consists of 
terms on the principal diagonal and the first term on 
each side of the principal diagonal. This solution uses 
minimum computer core. This type of matrix results, 
for example, from a one dimensional thermal model. 

Method 6 - This method uses a K-diagonal matrix which consists of 

the principal diagonal terms and K terms on each side of 
the principal diagonal terms. The value of K (the half band- 
width) may be specified by the user or automatically com- 
puted by FASTEMP. FASTEMP will also determine the minimum 
storage for the heat balance equations and set the solution 
method to 6, if advantageous. 

Method 8 - This method uses the Gauss-Siedel iterative solution 

The only matrix terms stored are the principal diagonal 
terms, nonzero off diagonal terms, and the constant terms. 

This solution will use a very small computer core space 
for many problems compared to the other methods discussed 
above. No computations are performed on zero off diagonal 
terms. The relative economy of using this solution will 
increase with problem size. The method is not recommended . 
for models with zero or low mass nodes. 
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4. PROGRA71 INPUT 

The input for a FASTEMP analysis consists of three distinct types; Job 
Control Language (JCL) cards, subroutine specification cards, and model input 
data. Figure 1 shows the order of the input card groups for a job that con- 
sists of one case. Data input forms are given in Appendices A and C. 

The JCL input is shown in Figure 1 as groups 1, 4 and 12. These cards are 
a requirement of the computer system being used. The JCL card input is not 
presented in this user's manual since the computer system requirements are 
frequently changed. The user should contact the computer support group for 
the current JCL input requirements. 

The user is allowed to code FORTRAN subroutines for inclusion in EQDAT calls. 
These are input as shown in Figure 1 as group 2. Group 3 is Subroutine EQDAT which 
calls the heat transfer subroutines, e.g., SRA (heat storage), DRA (conduction), etc. 

The model input data is shown in Figure 1 as groups 4 through 10. This 
data describes the model being analyzed and is explained in the remainder of 
this section. Group 11 is the last FASTEMP data card. 

4.1 PROGPvAM SPECIFICATION CARDS 

FASTEMP uses a dynamic loading technique which loads into the computer core 
only those subroutines required for each job. The Thermal Analyzer Subprogram 
is called by the main program and in turn requires u$er coding o^ subroutine EQDAXt 

The heat transfer method subroutines described in Appendix A must be 
called from Subroutine EQDAT. The subroutines supplied are independent and may 
be called in any order. However, the input data for the heat transfer method subrou- 
tines must be in the same order that the subroutines are called in Subroutine 
EQDAT. Failure to do this will cause a program error and a message to be 
printed. An example of Subroutine EQDAT is shown below: 


Card Column 

Contents 


7-22 

SUBROUTION^EQDAT 


7-15 

CALL*^QCRA \ 


7-14 

call’^sra / 

as many subroutine 
calls as required 

7-12 

RETURN ' 


7-9 

END 
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All calls to heat transfer method subroutines start in column 7 and 
are standard FORTRAN cards. Heat transfer method subroutines may be 
called more than once in Subroutine EQDAT but they must be separated by a 
call to a different subroutine to correctly separate data groups. Data input 
format options are grouped and may be a reason for calling a subroutine more 
than once as explained in Section 4.4. 

The card names must be punched in the card columns shown. Integer values 
on the case cards are indicated by a single card column number and must be 
right justified. Real values are designated with two card column numbers and 
must be punched with a decimal point and within the field shown. 

4.2 CASE DATA CARDS 

The CASEA card specifies a title for the case which is printed as the first 
line on each program output sheet. The CASEA card format is: 


Card Column 

Name 

Value 

1-5 

CN 

CASEA 

7-80 

TITLE 

User title 

The CASEB card contains 

the alphanumeric case name, the number of nodes 

print control, an orbit 

control 

option, the solution method number, and 

nonsequential node option. The 

CASEB data card format is: 

Card Column 

Name 

Value 

1-5 

CN 

CASEB 

10-15 

NCASE 

Alphanumeric case name 

20 

NCA 

Number of nodes 

25 

MPRNT 

Master print control 

30 

N^RBIT 

Number of orbits 

35 

IIETHOD 

Finite difference solution method 


3 - backward difference (square matrix) 

4 - backward difference (tridiagonal matrix) 

6 - backward difference (K-diagonal matrix) 

8 - backward difference (Gauss-Seidel 
iteration) 

50 IISvl Nonsequential node option 

0 - use sequential nodes 

1 - use nonsequential nodes 

The master print control may be used to surpress or to print subroutine 
and curve data input. fll'KfJT is used in conjunction with a print control on the 
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subroutine and curve data cards (NPRNT) . The value IIPRNT + NPRNT is computed 
for each subroutine and each curve and the following criteria is used: 
if (MPRNT + NPRNT) > 0 subroutine or curve data not printed 
if (MPRNT + NPRl^T) 0 subroutine or curve data printed 

The orbit control option, NORBIT, is used to recycle the program through the 
time history NORBIT times with initial temperatures continually updated. This is 
used primarily to simulate orbit operation but can also be used to calculate 
steady state conditions. 

Cases may be run with either sequential or nonsequential node numbering. 

The sequential option requires the nodes to be numbered 1 through NCA and the 
node numbers and equation numbers are the same. The nonsequential node option 
allows node numbers of 1 through 9999 and the equation numbers (1 through NCA) 
are assigned in the order that the nodes are encountered in the input data. 

The CASEC card contains the case beginning time, the case final time, the 
maximum temperature allowed, the minimum temperature allowed, the initial 
temperature option, a constant to be added to all initial temperatures and a 
solution method parameter. The CASEC data card format is: 


Card Column 

Name 

Value 

1-5 

CN- 

CASEC 

11-20 

BTIME 

Beginning time 

21-30 

FTIME 

Final time 


31-40 

TMA}{ 

Maximum temperature allowed (°R) - any 
computed temperature above this value will 
cause termination of the case 

41-50 

TMIN 

Minimum temperature allowed (°R) - any 
computed temperature below this value will 
cause termination of the case 

51-60 

TEQAL 

Initial temperature option if TEQAL 
> 0 - set all initial temperatures to TEQAL 
= 0 - read initial temperatures from CASET 
cards 

< 0 - use final temperatures from preceding 
case as initial temperatures 

61-70 

AODZl' 

Constant which is added to all initial 
temperatures 

10 
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CASET cards are only required if the value for TEOAL on the CASEC card is 
zero. An initial temperature must be input for each node, five temperatures 
per card. For a sequential node case, the temperatures must be input in the 
order 1 through NCA. For the nonsequential use, the initial temperatures must 
be input in the order of the equations. The equation order is a result of the 
order in which the program encounters the node numbers in the user's subroutine 
data and care must be exercised when modeling change cases which add or remove 
subroutine data cards and would change code number appearance order. 

A number of cases may be run serially with the final temperatures of each 
case used as the initial temperatures for the following case. The number of 
nodes in the following case must be less than or equal to the number of nodes 
in the preceding case. A larger number of nodes in the succeeding case will 
cause an error to be counted, values to be assumed for the missing temperatures 
and one time step to be run. Initial temperatures will be assigned using the 
equation order from the previous case. 

The CASET data card format is; 

Value 
CASET 


Initial temperature values 
(input in equation number order) 


Card Column Name 


1-5 

11-20 

21-30 

31-40 

41-50 

51-GO 


CN 


TI 


The case data is checked for errors, during the Z-phase, and errors may 
cause the following messages to be output; 

1) number of nodes wrong 

2) method wrong 

3) matrix half bandwidth wrong 
beginning or final time wrong 
T 


or T . wrong 
max min 

too few CASET cards 


4) 

5) 

G) 

7) too few temperatures from last case 

Missing or. incorrect data values are assumed and counted as errors. The assumed 
values which are used may cause or hide other errors. If errors are detected, 
the program will run only one time step. 
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4.3 DIITENSION TABLE CARDS 

The dimension table is used to input values which are required by the heat 
transfer subroutines. The dimension table input is optional since it need not 
be used for heat transfer subroutines data input. (See discussion of subroutine 


data input - Section 4 

.4.) The 

dimension table is input by a group of cards. 

The first card in the 

group must 

be a DIMENS card which has the following 

format : 




Card Column 

Name 


Value 

1-6 

CN 


DIMENS 

10 

NVAL 


Number of dimension table values 

The remaining cards in 

the group 

may be either VALUE or VALUES cards. The 

VALUE card specifies an index 

and one dimension table value. The VALUE data 

card format is: 




Card Column 

Name 


Value 

1-5 

CN 


VALUE 

10 

IND 


Index of Value 

11-25 . 

VAL 


Value 

The VALUES card specifies the 

value of four consecutive dimension table values 

The VALUES data card format is 

• 

• 


. -Card Column 

Name 


Value 

1-6 

CN 


VALUES 

10 

IND 


Index of first value 

11-25 

VALl 


Value IND 

26-40 

VAL 2 


Value IND+1 

41-55 

VAL 3 


Value IND+2 

56-70 

VAL 4 


Value IND+3 


Dimension table values may be input in any order. If a value is input 

more than once, the last value specified will be used for the case. If a table 

value index is out of the range specified for the table, that index 

will extend the table by as many values as required only when IND > NVAL from 

DIltENS card. Dimension table vaJ.ues which have not been defined will contain a 

322 

large number (1.265 x 10 ). Tlie card columns indicated for the number of 

dimension table values (NVAL) and the value index (IND) specify the right hand column 
of a four digit integer field, and digits in these fields must be right justified. 

The data card columns for the tal)le values represent real data fields and 
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values in these fields should be punched with a decimal point with right 
justification necessary for exponent (if used). 

4.4 HEAT TRANSFER METHOD SUBROUTINE CARDS 

The heat transfer method subroutines are used to describe the model being 
analyzed. Subroutines are available to compute various modes of heat transfer, 
e.g., heat storage, conduction, radiation, convection, heat flux, and fluid 
flow. These subroutines are available with rectangular, cylindrical, conical and 
spherical geometries to facilitate the description of the model. 

The subroutine data is input in groups of cards.. One group of subroutine 
data is input for each subroutine in the same order that the subroutines are 
called in Subroutine EQDAT. A group of data input for a heat transfer method 
subroutine must start with a zero card number. The zero data card format is; 


Card Column 

Name 

Value 




1-6 

SN 

Subroutine 

name 

(SRA, 

DRA, etc.) 

14 

NCARD 

0 




22 

NPRNT 

Subroutine 

input 

data 

print control if 


(NPRNT + MPRNT) 

■ ' T' ' > 0, input data not printed 

^0, input data printed 
26 HFMT Format if NFMT = 

^ general format 1 

2 general format 2 

30-70 IVAL Integer data for some subroutines - 

4 column fields 

The subroutine name must be punched starting in card column one. All other 
values on this card are integer values and must be right justified in the card 
columns sho\im. The remaining cards included in the group of data input for 
each subroutine are input with either General Formal 1 or General Format 2. 

General Formal 1 provides for the input of decimal data directly on the 
subroutine data cards. The data card format for General Format 1 is: 


Card Column 

Name 

Value 


1-6 

SN 

Subroutine name 

(SRA, DRA, etc. 

10 

CC 

Change code (See 

Section 4.9.3) 

14 

NCARD 

Card number 
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Card Column 

liame 

Value 

18 

NXl 


22 

NX 2 


26 

NX3 

Integer data fields 

30 

NX 4 


31-40 

XI 


41-50 

X2 


51-60 

X3 

Real data fields 

61-70 

X4 



The subroutine name must be punched starting in column one. Integer values 
are indicated by a single card column number and must be right justified. 

Real values are designated with two card column numbers and must be punched 
with a decimal point within the field shown. (Any exponent is right justified.) 

General Formal 2 requires the use of the dimension table described in 
Section 4.3. Only integer values are input on the card. Some of the integer 
values refer to locations in the dimension table where the decimal data is 
input. The data card format for General Formal 2 is: 


Card Column 

Name 

Value 

1-6 

SN 

Subroutine name (SRA, DRA, etc.) 

10 

CC 

Change code 

14 

NCARD 

Card number 

15-70 

NX1-NX14 

14 fields with 4 columns each 


The subroutine name must be punched starting in column one. All other values 
on this type of card are integers and must be right justified in their data 
fields. 

Card numbers are not required for a base case or change case; however, if 
individual cards are to be altered or added in a change case, card numbers 
must be used to identify the cards to be changed or added. No order of card 
numbers' is required and any number of cards may use the same number. All 
cards with the same number are dropped in a change case, if requested. 

All subroutines are available in General Format 2; however, the number of 
subroutines available in General Format 1 is considerably smaller. Some heat 
transfer method subroutines require more data input than is possible on a 
single card using Gener.i I Formal 1 . 
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The use of General Format 1 has the advantage of having the numeric values 
on the subroutine input card where they may be easily checked. An advantage in 
using General Format 2 is that the same dimension table location may be used for 
a series of node dimensions or multipliers and changing one value in the dimension 
table effectively changes these values used throughout the case. A model may 
be set up in such a way that both general formats are desired for the same 
subroutines. This requires that the subroutine be called twice in Subroutine 
EQDAT and be separated by another data subroutine. A complete description of 
the heat transfer method subroutines is given in Appendix A. 

The heat transfer method subroutine data is checked for errors, during the 
Z phase, and. the following messages may be output: 

1) Node number wrong 

2) Node connection wrong 

3) Subroutine references wrong dimension table value 

A) Subroutine called had no data. 

Additional information is output after the heat transfer method subroutine 
data has been processed in the Z phas.e: 

1) A nonsequential node case will output an equation number /node number 
table v/hich allows the equations to be identified. This is required 
to check out matrix error messages 

2) - A table listing the number of references for an equation number 

3) A table listing the number of internodal connections for each equation 

A) The total number of internodal connections 

5) The maximum internodal connection band (the matrix half bandwidth) 

6) A table listing the minimum equation number that is coupled to the 
equation listed 

7) A table listing the maximum equation number that is coupled to the 
equation listed. 

A. 5 CURVE CAPdJS 

The curves required for a case are each identified by a curve relative number 

(NRELN), a type (NTYPE) , and a class (NGLAS) . Curve relative number is generally 

assigned by the user, and the type and class are denoted according to the usage of 
the curve. Throughout this manual curves are designated by a three number system, 
i.e., "(NRELN, NTYPE, NCLAS)." The curves are a functional form which consists of 
one independent variable as a function of 0 through 2 independent variables (1 
through 3 dimensional curves) . The number of dimensions and the independent vari- 
ables for a curve are selected by the user. All independent variable values must 

15 

AffOOOMMfLf. oouat ..0 COHr*Of9ATIOM 



FASTEMP 


REPORT MDC A3791 


be computed before a curve may be used. Commonly used curve arguments are set up 
automatically in the program; however, the user may specify and set up special 
curve arguments. 

The curves required for a case may be input by two methods; same curve and 
temporary curve. Any required curve may be input by any method. 

The same curve method specifies that the required curve be the same as 


another curve. Input 

of this 

second curve by temporary curve method will 

satisfy the input requirement 

for the first curve. The SAMECR data card format 

is : 



Card Column 

Name 

Value 

1-6 

CN 

SAMECR 

10 

NRELN 

Required curve relative number 

14 

NTYPE 

Required curve type 

18 

NCLAS 

Required curve class 

30 

NREU'12 

Relative number of same curve 

34 

NTYPE 2 

Type of same curve 

38 

NCLAS 2 

Class of same curve 

The temporary curve method inputs the required curve directly by cards. 

The curve is input by 

a set of 

consecutive cards. The first card must be the 

curve identification card (TCIDNT). This card specifies the curve identifies- 

tion, the curve form and other 

information which describes the curve. The 

TCIDNT data card format is; 


Card Column 

Name 

Value 

1-6 

CN 

TCIDNT 

10 

NRELN 

Curve relative number 

14 

NTYPE 

Curve type 

13 

NCLAS 

Curve class 

22 

— 

— 

26 

NPRNT 

Print control if (NPRNT + MPPJiT) 



> 0, input data not printed 



< 0, input data printed 

30 

NDIMN 

Curve "dimension," may be values (0 to 3) 
(see discussion below) 

34 ) 



42 1 

NPTSi 

Number of points argument i (i - 1, 2) 

38 ) 



46 i 

NARGi 

Argument i type (i = 1, 2) 
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The second card is the curve title card (TCTITL) which is optional. The TCTITL 
data card format is: 

Card Column Name Value 

1-6 CN TCTITL 

7-80 CT Curve title 

The remaining cards are the curve value cards. The curve value data card format 

Name Value 

CN Blank field 

CV Curve values 

For curves with a value of NDDfN = 1 or 0, the dependent value is a 
constant, i.e., a function of no independent variable. The constant value is 
input in the second data field of the value card as shown in Figure 2k. 

Curves which are a function of one independent variable, y = f Cx) , have 
a value of NDIMN = 2. The NDIMN = 2 option is a curve which is interpolated 
linearly between points. The curve values are input with two point pairs per 
card as shown in Figure 2B. 

All curves with multiple independent variables must be input with a 
regular grid, e.g., an NDI1C7 = 3 curve is represented as y = fCx,z) and the sets 
of X values must be the same for all values of z. The three dimensional curve 
values are input in the following manner: 1) starting on the first curve value 

card, the values for independent variable one are input in ascending order, 
four per card; 2) starting on a new card, the values for independent variable 
two are input in ascending order four values per card; 3) starting on a new 
card, the values for the dependent variable are input in the order of all 
values for the second independent variable are input in the order of the first 
independent variable, all values for the second independent variable and the 
second value of the first independent variable, etc. An example of the input 
for a three dimensional curve is shown in Figure 2C. 

During the loading of a curve, extensive checks are made. Any errors 
which are detected are written on the output. These messages consist of the 
following: 
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1) A curve has an invalid identification 

2) A curve has been input which is not required 

3) A curve has been input more than once 

4) The same curve requested is missing 

5) A curve has an invalid dimension 

6) A curve has an incorrect number of values 

7) A curve has incorrect arguments 

3) A curve has independent variables which are not in ascending order 
9) A curve requested is missing. 

If any errors are detected, an error message is written specifying that the 
curve has been rejected. 

The card columns for the data fields of the SAMECR and TCIDNT cards 
specify the right hand card column of a four digit integer field. The data 
values must be right justified in these fields. The card columns for the curve 
value cards specify a real data field. These values should be punched with a 
decimal point with any exponent right justified. 

The user must specify the arguments used to perform a table lookup. The 
independent variable arguments are specified on the TCIDNT card. The following 
commonly used argument numbers are specified by the program: 

Argument Number Argument 

1 Time 

Method 3 through 8 - end of time step 

2 Temperature (can be the temperature of a node 

or other values - see heat transfer method 
subroutine writeups in Appendix A) . 

650+NN Temperature of Equation M. 

Other arguments are automatically computed for the heat transfer method subrou- 
tines and these arguments are defined in the subroutine writeups (Appendix A). 

The user may set up arguments for a specific use, e.g. , altitude, velocity, 

Mach number, etc., by using Subroutine Curves (see Appendix A). It is recom- 
mended that special - arguments be limited to argument numbers 11 through 19 to 
avoid conflict with built-in curve arguments. 

Two curves are always required for a case: 1) the calculation interval or 

time step curve (1, 10, 1), and 2) the print interval for outputting results 
(1, 10, 2). The input of these curves may be satisfied by either of the two 
methods of inputting a curve, i.e., temporary curves or same curves. 
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GENERAL HEAT TRANSFER 

TEMPORARY CURVE INPUT 2 

H 


MTTPE 

14 

NOAS 

18 

22 

Hi 

hCIVM 

30 



■i 




SFTS4 

58 
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7 

7 

7 



1 

HI 


IH 

■1 
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IH 




1 6 7 80 

TCTITL 

SAMPLe“cURVE°1D^ = CONSTANT 
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FIGURE 2A - SAMPLE CURVE INPUT FOR NDIMN = 1 


19 

MCDONNEL.L. OOUGt.A» COHr*OnATIOM 

















FASTEMP 


REPORT MDC A3791 


FASTEMP 

GENERAL HEAT TRANSFER 

TEMPORARY CURVE INPUT 2 



MTYPE 

14 

NO-AS 

18 

22 

la 

FOlhN 

30 

NPTSl 

34 




H 


NFTS 4 

58 

H 


Hi 


TCIDNT 

9 

9 

9 



2 

5 

■i 










ll 6 7 80 

tctitl 

3 AMPLE'^CURVEb 2 DbY = f(x) 

II 25 26 40 41 55 56 70 
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FIGURE 2B - SAMPLE CURVE INPUT FOR NDIMN = 2 
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FIGURE 2C SAMPLE CURVE INPUT FOR NDIMN = 3 
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4.6 CASE END CARD 

A case end card should be the last card in each case. The case end data 


card 

format is: 




Card Column 

Name 

Value 


1-8 

CN 

END^CASE 

4.7 

COMMENT CARDS 




Comraent cards may be placed anywhere in the user’s data for convenience in 
identifying data. These cards will not appear in the subroutine output and are 
only printed out by the input sorting section. The comment data card format is: 


Card Column 

Name 

Value 

1-6 

CN 

“h 

7-80 

CMMT 

User's comments 


4.8 MULTIPLE AI>1D CHANGE CASE CONTROL CARDS 

The program provides for multiple cases to be run in a single job sub- 
mittal. The first case run is assumed to be a BASIC case and no BASIC card is 
needed.' Parametric analyses may be run by inputting a BASIC case and then run- 
ning several change’ cases’ which vary desired parameters. Three types of cases 
may be run after the first base case, i.e., a BASIC case, a CHAI'JG^BASIC, and a 
CHAfJG^CHANG case. These cards immediately follow the END^CASE card of the 
preceding case. 

4.8.1 Basic Card - The BASIC case control card is used to input a full 
set of data which is independent from any previously run case. The BASIC data 
card format is; 

Card Column Name Value 

1-5 CN BASIC 

4.3.2 Change Basic Card - The CHAI'IG^BASIC case control card must be 
used to input changes to a BASIC case. The changes inserted for this change 
case always modify the last BASIC case which was input. The CHAi7G^BASIC data 
card format is: 

Card Column Name Value 

1-11 CN CIIAIJG^BASIC 

4.8.3 Change Cliange Card - Tlie CHANG^CHANG case control card must be 
used to input changes to the last change case which was input. A CHAI7 g'^CHANG 
case may be used after a CHAI>IG^BAS LC case or CHAI'Ig'^CHANG case. The CHANG*^ CHANG 
data card format is: 
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Card Column Name Value 

1-11 CN ceang’^chang 

4.9 CHANGE CASE DATA CARDS 

Change cases may input data to modify any of the case input in groups 6, 

7, 8, and 9 of Figure 1. Heat transfer method subroutines which were not called 
in Subroutine EQDAT cannot be added in a change case. The input data card 
groups for a change case must be input in the same order as the base case. To 
provide for the use of a subroutine which is not required in the base case, a 
call to that subroutine may be inserted in Subroutine EQDAT, and only a zero 
card input in the proper order in the base case. Data cards for that subrou- 
tine may then be added in the change case. The restriction of using only those 
subroutines called in Subroutine EQDAT applies to all cases which are input for 
a job submittal. 

4.9.1 Case Data Changes - Any of the case data cards may be altered in 

a change case by inputting the new case card in the change case. Any card not 
present in the change case data V7ill be taken from the base case. All CASET 
cards as well as a CASEC card must be input in a change case if any of the 
initial temperature values are to be changed by CASET cards. 

4.9.2 Dimension Table Changes - Any of the dimension table data may be 
changed. in a change case. The number of values in the dimension table may be 
increased or decreased in size in a change case. Dimension table values may be 
input xjith either the VALUES or VALUE cards. The change case cards are merged 
after the dimension table values from the base data, therefore values input in 
the change case will replace values from the base case. 

4.9.3 Heat Transfer Subroutine Changes - Any of the data for the heat 
transfer method subroutines may be changed. Temperature nodes may be added or 
deleted or the data for base nodes may be changed. The changes input for a 
change case are merged with the base case data using the change code in column 
10 and card numbers in column 14. The change code is ignored in a BASIC case. 
There are four change options available which have the follox^ing change codes. 

Change Code Function 

X The change code X ^^7ill delete all cards for a subrou- 

tine and must be input on a zero card in the change 
case input. 

R The change code R indicates that the card(s) from the 

base case will be replaced by the card input in the 
change case. For this option the card number must 
be nonzero. 
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Change Code Function 

D The change code D indicates that the card(s) from the 

base case with the proper card number will be deleted 
from the change case. 

"blank" The blank (and other characters not noted above) change 

code field will add the card from the change case 
input. 

Change cards for the heat transfer method subroutines must be input in the same 
general format as those in the base case unless a new zero card (using X for 
change code) is input to change the format. Change cards may be input in any 
order in subroutine data cards with the exception of the zero card, which must 
be the first card if it is changed. 

4.9.4 Curves Changes - Curves input as change case data may either be 
added to the change case or may replace curves in the base case. No provision 
is allowed for deleting a curve and it will be labeled as a curve which is 
not required, if it was replaced. A complete curve must be input for temporary 
curves, i.e., the TCIDNT, TCTITL and all value cards. No provision is made for 
changing individual vaj.ues in a curve. 

4.10 JOB END CARD 

The last card in a job submittal is the ENDJ0B card. This card informs 
the program that the last case has been read. The ENDJ0B data card format is: 

Card Column Name Value 

1-6 CN ENDJ0B 
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5. PROGRAM OUTPUT 

The output of the Thermal Analyzer Subroutine consists of three types of 
information: 1) output of the input data; 2) output of the results; and 3) output 

of error and warning messages. Most of the output of the input data is optional 
as specified by the user. The method of controlling the output of input data is 
discussed in Section 4. Output of the program contains some sections which are 
always printed, however , the majority of the output of the results of a case is 
under the user's control (e.g. node temperatures may be output at any problem time). 
Output of node temperatures is controlled by the print interval curve (1, 10, 2). 
Error messages are output specifying incorrect input data. Errors found in the 
input data cause the program to compute one time step only and to terminate the 
case. (See Section 7 for an explanation of the procedure used.) 

The first page of output gives the computer time at the beginning of the 
M-phase and the number of core locations available for storing the data for the 
case and the matrix. 

The next page of output for each case gives the computer time at the start 
of the case, the date, and the subroutines called by Subroutines EQDAT and S0LVE 
which require input data. Cards which have been rejected in the M-phase processing 
are also printed on this page (see Section 2). This output is continued on 
succeeding pages as required. 

The next page of output gives the case specification data from the CASEB and 
CASEC data cards. The title from the CASEA card and the case name from the CASEB 
card replace the standard page header title and case name on this page and all 
following pages of output. Any errors in the case data will cause messages to be 
printed on this page of output. 

The next page of output gives the dimension table if one is used. Any unde- 
fined dimension table values will contain the number 1.265 x 10^^^. 

Tlie next output from the program is the input data for the heat transfer 
method subroutines called by Subroutines EQDAT and S0LVE. This output is optional 
and the user may specify individually which subroutine data will be printed. The 
output for each subroutine starts on a new page and continues on subsequent pages 
as required. Errors in the heat transfer method subroutine input data will cause 
messages to be printed on the output immediately after the card containing the 
error is processed. Any CDC hardware error codes occur for the next card which 
is not printed prior to error and should be noted when debugging jobs. 
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Following the processing of the heat transfer subroutine data in the Z-phase, 
several tables are output which contain information about the heat balance matrix. 

If a nonsequential node number case is input, a table is output of node numbers 
versus equation number. Tables are output which give: the number of times each 

node number is referenced by a heat transfer method subroutine; the number of 
internodal connections for each node; and the total number of internodal connections. 
The maximum internodal connection band or matrix half bandwidth is also output. 

Two tables are printed which contain the minimum and maximum equation that a given 
equation has connections in the thermal model. 

The next section of the output is the curve input data. This output is also 
optional and the curve data printed may be specified by the user. Same .curve 
messages are output with the temporary curve data. Error messages will be output 
for curves which are incorrect and have been rejected. 

After the curve data output, two messages are output giving the number of 
locations available before and after the core was requested for the matrix. The 
time used to the beginning of the P-phase is also output. 

The results of the problem solution may be output as specified by the user. 

The temperatures for each node may be output at any problem time. Some heat trans- 
fer subroutines have additional output which, if desired, is printed at each 
computed time step. Tables of maximum and minimum temperatures for each node and 
the problem times for which they were computed are output after the last time 
step is computed. 

The next page of output gives the P-phase computer time elapsed for the case, 
the number of time steps computed and the number of errors detected. This is 
normally the last page of output for a case. 
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6. SAMPLE PROBLEM 

The sample problem presented in this section illustrates the setup of a 
model for a typical heat transfer problem. The sample problem does not contain 
all of the possible combinations of program options which are available, however, 
enough of the options are used to allow the user to gain a familiarity with the 
general setup of the deck. 

The sample problem is a three dimensional heat transfer analysis which has 
both steady state and transient computations. Modes of heat transfer considered 
include heat storage, conduction, radiation, convection, and interface 
conductance. 

The model used in this analysis is shown on Figure 3. It consists of a 
symmetrical section of the lower insulated surface of an actively cooled 
hypersonic aircraft external skin panel. Aerodynamic heat transfer through the 
insulated surface is absorbed by methanol /water coolant circulated through Dee-shaped 
tubes bonded to the inboard surface of the skin. Coolant temperature at the 
panel inlet was assumed to be 0°F and coolant flowrate equal to 77 Ib/hr per 
tube. The panel skin and coolant tube were assumed to be constructed of 

3 

aluminum with the back-up material of 4 Ib/ft Johns-Manville Microfiber insula- 
tion, and the outer skin of titanium. The basic model dimensions (i.e. panel 
length, tube diameter, tube pitch, and material thicknesses) are shown in 
Figure 3, with corresponding node number identifications. A total of 208 nodes 
was used. The nonsequential node numbering option was used, which allows 
identification of different materials or sections of the model with groups of 
node numbers. 

All nodes of the model are three-dimensional and are described by rectangular 
geometry, with the exception of those for the coolant within the tube, which are 
defined by cylindrical geometry. Nodal dimensions as well as values required by 
the heat transfer subroutines v/ere input on the dimension table. Figure 4 
illustrates the dimension table input parameters for the sample problem. The 
variable terms shown on the dimension table refer to model node dimensions and 
coolant Reynolds number transition from laminar (Re ) to turbulent (Re ) flow. 

Use of the dimension table, as illustrated here, provides a considerable degree 
of flexibility for implementing model or subroutine changes. 
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The mission profile used for this sample problem included an initial steady 
state cruise condition at Mach A. 5 and an altitude of 90,000 ft. During this time 
flowrate through the coolant tube was maintained at 77 Ibm/hr, with a panel 
inlet coolant temperature of 0°F. At a time of -15 seconds, a total failure of 
the cooling system was imposed. This was followed at a time of 0 seconds by 
initiation of an abort descent trajectory resulting in lower aerodynamic heating 
rates. It was assumed that the cooling system failure was of a nature which 
caused instantaneous depletion of the coolant in the tube. Aerodynamic heat 
transfer rates during cruise and the abort trajectory were input in terms of 
external surface heat transfer coefficients and adiabatic wall temperatures 
versus time. These were calculated separately from other analyses based on 
aircraft geometry and parameters of flight altitude, velocity, and angle of 
attack versus time. Two cases were run for the sample problem. The first case, 
or base case, calculated coolant and panel structure temperature distributions 
and coolant pressure drop during cruise prior to cooling system failure. The 
second case, or change case, calculated panel transient temperatures versus 
time from cooling system failure (-15 seconds) through the abort trajectory 
(600 seconds) . 

Initial temperatures for the base case were input for each node as an 
example of the CASET card input. This was not necessary for this problem since 
steady state equilibrium temperatures were computed by starting the case at a 
time of -300 seconds and computing time steps with no heat storage until time 
equals -15 seconds. (See Section 7 for a discussion of the method used to 
compute steady state temperatures.) Temperatures were computed at the calcula- 
tion intervals of 30 seconds for the base case and 5 seconds for the change 
case as given in curve (1, 10, 1). Print output of all node temperatures is 
controlled by curve (1, 10, 2) and was performed at 100 second and 15 second 
intervals for the base case and change case, respectively. Material property 
values (e.g. density, specific heat, thermal conductivity, etc.) as well as 
interface conductance values were input in curve form either as constants or 
versus time or temperature. Curve inputs were also used to describe other 
parameters required by the heat transfer subroutines such as external surface 
heat transfer coefficients, adiabatic wall temperatures, coolant flowrate, 
coolant friction factor, etc. A total of 29 curve inputs was used for the 
sample problem. 
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Heat storage was computed for all nodes except those of the coolant by 
Subroutine SRA. Conduction was computed by Subroutines DRA, DRB, and DRC. 

DRA was used to compute a conduction heat flux between nodes of the same 
material (e.g. nodes 1 and 11). DRB was used to compute a conduction heat 
flux between two nodes of different materials which are separated by an inter- 
face conductance (e.g. nodes 1 and 301). Subroutine DRC was used to compute 
conduction between two nodes of the same material which are separated by an 
interface conductance (e.g. nodes 1 and 101). Aerodynamic heating to the 
external surface (nodes 601-639) was computed by Subroutine VRA, with heat 
transfer coefficient and adiabatic wall temperature inputs. External radiation 
heat fluxes from nodes 601-639 to a 59°F sink temperature were computed by 
Subroutine RRA. The forced convection Subroutine VCA was used to maintain 
the coolant at the panel inlet (node 200) at 0°F by inputting a very high heat 
transfer coefficient and a 0°F adiabatic wall temperature. Subroutine FLUID B 
was used for the base case to compute convection heat transfer between the 
coolant and the tube wall as well as pressure drop down the length of the tube. 
Laminar flow was assumed for this case to apply below coolant Reynolds numbers 
of 2100, with a step change to turbulent flow at higher Reynolds numbers. 

Figure 5 shows a listing of the input data cards. This listing includes 
the program specification cards (i.e.. Subroutine EQDAT) . The job control 
language (JCL) cards are not listed in Figure 5. Figure 6 shows the output . 
for the sample problem. Sheets 1 through 23 show all of the output possible 
from the Thermal Analyzer Subprogram for this sample problem during the M and 
Z phases. The user may partially or totally prevent the output of the sub- 
routine data (Sheets 5 through 15) and the curve data (Sheets 17 through, 21). 
The remainder of the output (Sheets 24 through 28) show the output from the P 
phase. A change case was setup for this model as shown in Figure 5, to run the 
transient analysis. Figure 5 shows the method to eliminate a complete 
subroutine's data (X in column 10 of FLUID B card). This simulates loss of 
coolant fluid at -15 seconds. Only part of the output for this transient 
is shown . 
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MCDOMNeui. oouoLAm cort^ortJkTiOM 


FIGURE 4 - SAMPLE PROBLEM INPUT DATA CARD LISTING (SHEET 
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MCDOMMmt,i. oouaLAm cortr^oftMnoM 


REPORT MOCA3791 


ORIGif^AL PAGE IS 
OF POOR QUALITY 


FIGURE b - SAMPLE PROULEM OUTPUT (SHEET 1) 
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AtCOOWWffLI. OOL/a<- A V COW^OMATtOM 


ORlGJf^AL PAGE !S 
OF POOR QUALITY 
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REPORT MDC A3791 
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mtcooMMmLL. oociOLAS com^ort/knoN 


onmukl PAGE IS 
OF POOR QUALITY 


FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 3) 
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JMCOOMMf LI. OOUOLJkU COttfOftMTiOM 


FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 4) 
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MOOOMiVfI.1. OOC/OLA9 OOMr#»OM^riOW 


0 RiGi? 4 AL PAGE [3 

OF POOR QUALITY 
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mCDOMMeLL. OOUGLAm COftr*OfWMTI€>M 


FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 6) 
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MC00MNK1.L. oouaLAm cowp»o«i * rioM 


ii#^ODWCffiILITy OF THE 
0^«:aiML PAGE IS poor 


FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 7) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 9) 



FASTEMP 

) 


) 


REPORT MDC A3791 


••*•••»■« * •••** »•»«>.• 'T'rr*. ;r^ — ■'•" u •••** •?• *' ~ '•'- . .. - . . .'v'j - ric •• !.•-«« « i t $ 

r »'■ ^•'Q »•■■»* rv«»'>i r f *••*'/■ f- c. j " ““'•‘m a r. "“«r f'.fvfv. Jii u'ltUiHt' 

CNjfVir «■ O O^O'O' OI.'IN. <\.(V U'** J J 0*i» tP«M ^ J J OO<SlfSj<\jU'«.'0'U'0'0’«;\£*il'**»^ 

.^^fvjfv€VjOw»r**nfoo<v»-J •# J j j M«i 4 r\,iNif\iC 3 ^p"'>p>orvjCT' 0 '( 7'«4 

• •••••••«••■•••••••••••••••••••••••••••••••••••••«•• •r?fwfk»P>tA 


It It II II M II II II II II II II II II li II II II M II HJ |t II II II II M II M II ft II It II l« II M II II II II II II II II II II II II II II II II ,1 II II It II 

>O>X>*>»X><’X>«>cXX)*^>'>«>>M>«M>>XMMtaX^>)kX>(k'X»XX>(X>(>i*(KXHXXMXXX> 0 «VHXXX 

C‘C-'c;c_.u:ccac-CiL-t::c.c:oc.-cco::rcr'CCiO»-coc-c cacc-ocuc ccjoccrcoccr.oc^cc c^c-ooooo 

'S.XXXXX'S.XXXXXV'S.'vXXVX 'S.'X'VXXXX.X XXX V>* X XX XXXXXXXXXXXXXXXXXXXXXXX 
u>t< I. 'UlU 'iilU -U -u > 1 . -l< 'L 't-. 'li -t. •litL 'UlU'UUiiHJ-U>UJUMtUi'UIU*lt OUl lU UHl^U 'U Li 'U 'U <li •U.UUilb 'U 'U.'tilU'U 'lilU )b 'U'itIliJU Hi I 

onaaaaaooooQO'OCkOnaoaiacfcQftan'auoonouooQonaaQficr&QOQocxouQorLoaaa: 

<t< 3 d<l'T«l« 3 <X<X« 3 <<l« 3 'X« 3 <X< 3 <X 4 <X<«<l«t< 3 <K 4 <3 <3 44 <<* 4 «l<l 44 <<I«l««a«i«X 44 c 3 «x «34444444 


oocacsooeo 0000000^0000 ooc*c 30 <^oooooooooc-eoocdooonooe>oooc 90 oooooo 
oooooooooooooooooooooeooooeoooooooooooooooooooooooeeoooooo 


w^oooooooC'ooeoooaoooeoec'ooooeCt.'oeeooooeooooeectoooooooooDoOoo 


x-oo ooonooeooooooooctooooooooooooOoooooooeooooooooocoBooooooo 

I I I « I « I I I I I I . I I I « I I I I I I • I I I 1 I I • I I • I I * I « I I I • I I I » I I I I I I I I I 1 I I 


«£ oooo OOOOOOOOOOOOOOOOOOCaOOOOOOOOOOCJOOOC OOOOOOOOOOOOOOOCJOOOO 


f^rsooooooooooooooooooooooooooooooooooeooeoooooooooooDooooooo 


■eOOOOOOOOOOOOOOOOOOOOOOOC OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
I 


.^ooouooooooeoooooooacioooooooooooeooeoanoccoocdoooooooooooooo 


O >0 000000000 00000000 0000000000000000 0000000000000001^0000000 

' o 


fV'C'rx'^a 7 ®'«xfo< 0 'C»-x»o-r:'r 5 -r;x ■oo'r>«x*no<C'*>fx*oo'C t'X i'o o cfx'o«-^T*x>o oi©crx^®'C‘r>x»<-'«o« ox 

»o*o»o»o' t» 7 toiii)i«> 


to 






5 


o 


OOO 0 cj 000 oa 00 0 00 0 ‘^i- 5 ' 0000000000000000 c* 0 o 0 f 3 c 30000 OOOOin OOOOOOOOCJO 




^t'n f*rfX-» tr ' 3 jC' 0 XT<T'^-^*'j"njt/''/'x«{yn t >f' i'Xo T'c>— 

I f . .-iT y |■'•J■^ XX X X •w. ^ f *r «ry*-rx kx x r rs,rv"vj<‘>t'vsr y-rc xxi-«r«CT"'j?‘vP.'p''^xxxN x<Qo«r,-«>rT« 


■n T «T' •» .f' Xxoj<T^ — 'j*o ^tTv.rx — <*j*o j •."> px nO'Cj-^rvj*^ ^ P\rx ■© <r- ca *o y lO^rx oP'o~4»vjf^ ^ip 


1 y •!» ^OXOiPO ^ir 0>x xa'0«-»l\f-0 CX !©•»' -j»-«l\|»0 £|X «0 J'-3 lO«£X O '}><-> ^INi»0 ^Ift.^x «<1* 

t * J ty T ^ r J t P P r P't* .tvp 'O PXXXXXXXXXX-O c 0'P"053T>TtO^ 

CJ 


•o 


I <> 'f 14.' *> < 1 . .'i cc (1 iir ^ a u .(.'.r a. a 3. it'u' o u. lira ^ a-x a o' x.o v a a: ora a u Kt-x a..x a a y a 3 . a a o. ixu 'x'oraa or or or 
-I'^ooTTDoaooo'Tinc'm'aoo “>ooo ‘0 -jooooto orir^fDO ;~-i.-p' 3 ^o' 3 n 3 n'Oo^)'T:i pjorsnooo 


2D 

O 

02 

Qu 


00 

I 

LiO 



58 

mucDOfy/vEi. l. ooocxxs contort myioiv 



MtitHAN'lNS HT Sm-lLj maCii -.K tN JM IN^.UL .25 tN. Tunr •• 75 / 12 / 0 ».» 20 . 52 . 01 . » P 45 E 9 - CAS-i MP 7 S £0 = 


FASTEMP 


REPORT MDC A3791 




lr>lr'u^u^u^u^u^u'a'u^tJ«lX^if^u^U'U'UkU«U'U'oc>r'«cJC^C'c^c- u>oc>oi.t« >(«(.^r>c'c:voc:ir'>oc>o«i.p«A>(>oooo 
C'*'-*0C.'00C'0000 tJ »JC* r_'f^CJO«aOC»Crc ^**-4 'C' <, 't^C. C.‘C. C'lA^U-'OO 

■ • • I •* •* • • •« ; • i» •» « I «rs' ~*i i 

t.’i.'i.>Ai 4 'i''U''**u A.i.-’u-M. L li uL-u ■ f ‘t-t 'e Q-» 4 u*o«o«u»»fv 

,jip>..ps.^>.ci 0 ’C' 0 'C*'t-p'^*> ® «•«*»>.. C'OC'C'C'C^OC 

LP po>i>vL'\ 4 '^^u>u' u't/'U* i/iiociu >u'u«u'y'u>u*^i.'* i.*c«rAc iocorr oo cjcjcjc^coc.* o«>< v 4 
lPr'(— >r'P' ■ ■ . . - - — 


f«~ p>. !>.. rk. ^ ^ ^ ^ w-i «-• 


Qi C'O'AOO. *r « 


II II It II II M II M II II (I II II M If II II It If It II It II II II If II II II fl K II « If II II f« fl II II H II II II If II ri If II II II II If II fl II n H 
XXX»«ik'>«X> 4 X>X»X>«>«>«X>«>i>«>«»*VV^X>'>l>«'*»>»X>t>.v> 0 «XX»r>»>l>*>*>i>i>»>>rXXX>eXX»«>«>r 

c?cc:ot::cjcc.c-oc:cocjcciL:ac;r5c.cit r.t-f >oc*c'c.;ccoc oc oc*ooc:t;GC OGGi_*c c:aGOC?ooo 

XV's,>»'s.V'XV>..XVv> 4 '«.V.'>..>..'S.»s.’vN.’s'*.*x>.'N.'» 4 's.N 4 ’V> 4 >»'.fcXV>*,V,XV*VVX«V>.'«>.V.X.’>SVX’**.VX.-V> 4 'VV*\ 
ItIU.Il'UlU UIU>U'U-ii>U*U.| 4 >li' 0 ‘(tlUitt'U*k>*U.U l••k>‘U UiU*UII.'U-il'U'i''lA;ti'lUU'll'li'U‘LUl'U<iJIUfUl 

c^ooooooi>oQounoaii.QOCkCOoc^uAAf>ooo 3 ^Qft'iyA'np(Nuooniiiifvaa.ooACk{ynoftUQar 


or«c 3 aooaaeoe 30 oc^r>c 9 oooi 3 c:Oaoc»ooooeoc>ooooooc«ooooODDooor>oc.>ooo^o» 

Illllll•l•llllllllllllllttttllllllllltll»ltilitttl|l|l|||| 

oejooooooociooaouooooooooooooooooooooooooc^eoDOOOooctoooowoooc* 

o 

« 


»<oe.,c>e^oo c 3 0 C*i^ot:* ooooooo*->C'«^oC*c'«'*ooo«oc»o» >000000 013 oO ffoe'ooooc'oociooo 
'^li•••«ll«llll«lll«llll«l•lll•lllllli•illllliltlllltlilllll 

o 

K. ooofjoi ic'ooo0cjoooooeior»i^oc»0f»ooc»r*ooc^or>ot>c>f »ooOc*oo Of *00^^00 t»oo 00 ou 


^.'coc oooaoeooeoooooooooooooooooooooooooooooeoo ocd&>ooooiAOt_>oo eo 

*0 I I I I I I I « I I I I I I I I I t I I I I I I I I I I I • I I I « I t 1 I I i I ( I I I I I > I I I I I I I I I I 
t 


Ajooooooooooooooooooooooooaooooooooooooooooooooooutoooooooo 00 
•C I I I I I I I I I I « t I t I I I I • I I I < » « « t I I • I I * I I « « > I I I I I I i I I I I I I I I I ) I I I 
I 


«cooOoooeoaoeoooooo« 9 oooocjooooooooeoooooeoooooooooouGooe»ooooo 

I 


■4 oanoooooeeooooeooooooaooaoooooooooooooooeoooooeneoooooooao 

I I I I « I I I t I I I I I I I I I I «« * I I I •• I I I I I • I I I I I t I I I I I I I I I ( I I I I I I I I I 
I 


osooooooooooooocoooooooooooooooooooooooooe OOOOOOOOOOOCOOOOO 
I 


• ^00 o 00 000a 000000000000 00 ooooooooooo 0000000 ooooo 00 

•4 *4 T" 1 PO — -n rO»n ^ j j ^ -J 


pg 0 - o o o o K r^> K. K. r». r>> rw fs. X. r>> ho N. K. K fw N. ^ ^ ^ ^ ^ ^ ^ ^ ^ >-4 rvj CM ^ 

■4 rrnm'n pt "opopo 


«" 'Vhf'M (NjAjhJfyiMfMINjtT'O^r'O ft" 0 \ 0 »r^ *• «r«n <'0 **«■*«-• o^»VIVIVlf<Jf\JfVfM(Vfv;PvjfMfV'M CM f'JfMfyiMf'jtf'UMfM/'UN 


t IT ir>:iM.‘p.kr (\iiVMN«rvi<M-^ ^ ^ o«'^^ir‘p^io*of>?*p>ro^pp>pootf‘«n cic <r:V«i«.a 4 ccw t, x: ^ 

o ,<» 4 « 4 o.i»-t»<wwo^p^'^p'^ro»opp>p>^."opo»p 5 *rp'>po*rr»*» 5 »p 5 pp 5 x>*o»»>’p>|p> 


>0000000000 o- >ci'nof 30 oo« Jf 5 oc>i.»c»ooooo .TirT*:>oooo 000 oooo 00000 0000 eSPs 

*~i>iiltiiiiltltitliii*iiili|i«il**i*iiililiiiiiliitiiliitli 

4 .T -I 4 4it fV«jfM<Vhj.*Mr\jrvjp»ji\|IM<M<MlMIM«MfSlhj»MIVCSjCVIMMIM 3 4 3 4 4 4 4 4 t -4 4 4 4 4 3 4 4 444444 


hi ^“>J*** f ir..Oh> «CT‘-»r.»*h O fh. flOCT' •^IMO 3 iT .Oh*p>P J •-•P\>ph msTh*. JT J'h«..r x-'^**i4 T I> -3 Tw ® 

fT*r»Q,p,c 

x.Ch»'nfT'o-^.*j'hj r> — >%(»•. jj> Oh-“n<»— 4 ivi*o^u' j'K. .<r-jp»* $ o'vf'p* -o ■hu*' ^ »jpsi*op>*. -c ccnjon- tj rsi **> - o ovT* 

c -r a- .0 r*n »npo"i**^t: r. t »> t. )'c» 4 i»»'.» o- »• x ^ §• .r* ^ ."M nj rv pj 'O -> *r» j j r r n'-h 

r c 5 *^hj*«o !ii..th. * 0 ' 0 '.'>o .to j> o i>ou>«-»'Ac>'.i>o /'oino«»tr'y>.ou'o A Oh. ^’ro-^fVPh^iAXh- 

c P>- ^ ■ >o^^P^*\tr;«MAjXi '•>«P»Pn*hi<^~t«0'p'OI'»J r J 3 4 3 -4 .3 

I hjr-j(vhirvj(sj-Mrghjrgh)inro>p)i«thm<p|*»o*p>'OPp>OPr*p)poo*toX t x TX 44 34 344444444 

u 


.0 


f -t-x f p-r 'p *x <t 'f ^ <t *1 <1 «T -r^r '1 'Y't't'i'f <f pi'f 1444 »t 4 i ♦^-i 4 4^7 4 <t4 «t4 

10.3: ku-i'.VixaxaxOrJia.Jtxoik.-x.Ttcrwa^ tic 3^ 0*0 u a ik.i u:a o.k u.> My a a a a. y > a o-x t.4 uiuocol 
*Do anniT"j'^oo^r 3 Tno 300 " 3 orv- 5 '>C“^no'^ 0 'T *'3 jr>o a- 30 *ooooii 


59 


MCOOMMBLt. OOUaLAm COttfOtWJCTIOM 


FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 11) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 12) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 13) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 14) 
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FIGURE 5 ^ SAMPLE PROBLEM OUTPUT (SHEET 15) 
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S?.^?^ODUCIBILJTy OF THE 
■ "mAL PAGE IS POOR 


FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 16) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 17) 
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oociOLAs comr*oftM tiom 


figure 5 - SAMPLE PROBLEM OUTPUT (SHEET 18) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 19) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 20) ' 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 25) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 26) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 30 ) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 31) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 32) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 33) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 34) 
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FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 35) 



FIGURE 5 - SAMPLE PROBLEM OUTPUT (SHEET 36) 
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7. PROGRAM USAGE 

Section 4 has specified the data input for a case by discussing the several 
types of data required. This section includes several additional requirements and 
general suggestions for usage of the program. The section on data input also 
discusses errors in the input data as it is sorted and read. This section discusses 
additional errors which may occur and what the computer program does in response 
to these errors. 

7.1 MISCELLANEOUS REQUIREMENTS FOR PROGRAM USE 

The (I, 10, 1) curve (calculation interval) and the (1, 10, 2) curve (output 
interval) are required for every case. An error is counted if these curves are 
omitted from the data input. (See paragraph 7.2 for an explanation of the proce- 
dure the computer program takes when a case contains errors.) 

Subroutine CURVES is used to perform table lookups on curves not specified 
on heat transfer method subroutine data input cards and to set up arguments which 
may be used for table lookups in curves specified in the heat transfer method 
subroutines. An example of this application is the altitude versus time and 
velocity versus time curves which are frequently used to analyze an aircraft 
mission. Neither of these curves are called for by the heat transfer method sub- 
routines directly, however, the information may be needed for the computations. 
Subroutine CURVES must be called in Subroutine EQDAT prior to calling any 
subroutines which may need arguments setup. Failure to call Subroutine CURVES first 
will cause a value of zero to be used on the first time step and may cause a 
division by zero error. (All undefined argument values are zero). The curve 
arguments may also be setup one time step late for the remainder of the case 
depending on the argument used. 

7.2 SOLUTION ACCURACY 

Solution accuracy may be influenced by several factors which are under control 
of the user. The backward finite difference solution does not alter the time step 
input by the user. Increasing the time step tends to Increase truncation errors 
for transient heat transfer problems. The higher the rate of change of tempera- 
ture with time, the larger the error is in the computed temperature at the end of 
a time step. One of the causes for this is the evaluation of properties at the 
last known temperatures. The linearization used for T^ in the internal radiation 
heat transfer subroutines also causes errors in the computed heat flux which 
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become larger as the rate of change of temperature increases. The user should 
specify the time steps for a case to reflect the heating rates input for the model 
boundary conditions. 

7.3 STEADY STATE SOLUTIONS 

Cases may be set up for steady state only or a transient analysis with steady 
state temperatures, optionally computed at one or more problem times. Cases of the 
first type may be computed by omitting all heat storage subroutines and running 
the case for several computations (time steps). This is necessary to allow 
temperature dependent properties such as specified heat and conductivity, to be 
evaluated at the steady state temperatures since they are evaluated at the old 
known temperature. Some heat transfer methods (e.g, internal radiation) employ 
an approximation used to linearize the heat balances. These subroutines may 
require a few additional iterations (time steps) to reach a steady state solution. 
Ten time steps will produce a steady state solution in most cases with no heat 
storage. 

Cases which combine steady state and transient portions obviously may not use 
the method discussed above. Heat storage terms must be included for the transient 
portion of the analysis and zeroed out for the steady state computations. This 
may be done by inputting a density curve as a function of time and specifying the 
density equal to zero for the time intervals when steady state temperatures are 
desired. The number of time steps required to reach steady state temperatures 
will vary as described above (size of step is immaterial). 

7. A USER ERRORS 

There are two categories of user errors; those which occur before the P-phase 
and those which occur during the P-phase. Errors in the input data are detected 
and error messages output in the M and Z-phases of the program. The computer 
program continues to process cards, regardless of errors which are detected, until 
the completion of the Z-phase. The program will compute one time step in the 
P-phase if errors have been detected in previous data processing phases. This 
is done to check as much of the data input as possible and therefore minimize the 
number of runs which must be made to obtain a solution. There are situations 
where one error causes a number of other error messages to be output. 

The computer program assumes values for missing or bad input data to allow one 
time step to be computed. A value of 1.0 is assumed for a specified curve that 
is missing. A value of 1.265 x lO^^^ is assumed for dimension table values that 
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have not been input. The program assumes the following values if the CASES card 

is not found t 

NCA = 10 
MPRNT = 0 

NORBIT = 0 

METHOD = 3 

The program assumes the following values for the CASEC variables if the card is 

not found: 

BTIME = 0. 

FTIME = 2.0 
TMAX = 1000 
TMIN =2.0 
TEQAL = 560. 

ADD2T = 0. 

Additional case data values are assumed if the following errors are detected: 

(1) The value of NCA must be greater than one for all methods. 

(2) The method number must be valid; otherwise it is set to Method 3. 

(3) FTIME must be greater than BTIME; otherwise they are reversed. 

(4) TMAX must be greater than TMIN; otherwise they are reversed. 

(5) TEQAL has three options. TEQAL = 0. causes initial temperatures to be 
read from CASET cards. Either too few or too many initial temperature 
values being input will cause an error. The additional values of 
initial temperature are set equal to 560°R when too few initial 
temperatures are input. If TEQAL is negative, final temperatures from 
the preceding case (if present) are used for initial temperatures 

for the present case. A value of 560°R is assumed for initial 
temperatures if there are less nodes in the previous case than in 
the present case. If TEQAL is greater than zero, no error checking 
is done. 

The following user errors occur in the P-phase. 

(1) Temperatures greater than TMAX or less than TMIN cause an error to 
be counted and the case is terminated. 

(2) Divide checks occur when the computations result in a division by 
zero and may cause mode termination by the CDC computer. Input data 
should be checked for zero values. 
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A warning message is output each time a curve is extrapolated. A curve 
argument outside of the range of the input values for the independent variable 
causes a linear extrapolation using two points closest to the value desired. 

This message does not cause an error to be counted as the value may or may not 
represent an error. The user should check the extrapolated values to determine 
their suitability. 
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APPENDIX A 

SUBROUTINE DESCRIPTION AND DATA SHEETS 

Section Page 

A. 1.0 Heat Storage Subroutines A-3 

A. 2.0 Heat Conduction subroutines A-14 

A. 3.0 Heat Convection Subroutines A-43 

A. 4.0 Heat Radiation Subroutines A-51 

A. 5.0 Heat Flux Subroutines A-72 

A. 6.0 Fluid Flow Subroutines A-80 

A. 7.0 Model Subroutine A-90 

A. 8.0 Miscellaneous Subroutines A-95 
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Curve Specification 

Almost all FASTEMP subroutines require as part of their input, curves which 
describe material properties, source temperatures, coefficients, and other boundary 
conditions. Section 4.5 describes the input and use of these curves. Curves are 
designated using a three number system (NKELN, NTYPE, NCLAS) . The following is a 
standard list of curve NTYPE and NCLAS. 


NTYPE 

1 

2 

3 

4 

5 

10 

10 


11 


11 

12 

12 

12 

12 

13 

13 

13 

15 

17 

19 

99 

99 


LIST OF CURVE NTYPE & NCLAS 
NCLAS CURVE DESCRIPTION 

0 Density (Ib/ft^) 

0 Specific heat (BTU/lb°F) 

0 Conductivity (BTU/hr-f t-°F) 

0 Emissivity 

2 

0 Interface conductance (BTU/hr-ft -°F) 

TIME STEP CURVES 

1 Delta Tau vs. Time (Seconds), NRELN = 1 

2 Temperature printout, NRELN = 1 
ATMOSPHERIC PROPERTIES 

1 Temperature (®R) 

2 Pressure (psia) 

TRAJECTORY DATA 

1 Altitude (ft) 

2 Velocity (ft/sec) 

3 Mach number 

4 Angle of attack (degrees) 

TEMPERATURE CURVES 

1 TAW - Adiabiatic wall temperature (°R) 

2 TR - Radiation Source temperature (°R) 

4 TP - Known temperature for conduction (°R) 

MISCELLANEOUS CURVES 

0 h-convective film coefficients (BTU/hr-f t-°F) 

0 Shape factor curves 

1 Fluid flow rates (Ib/hr) 

0 Heating Rate curves (BTU/lir- 1 L ^) 

0 Multiplier Curve 
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A. 1.0 Heat Storage Subroutines 


Section 
A. 1.1 


SRA, SCA, SKA, SKB, SKC, SSA 


Page 

A-4 
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A. 1.1 SRA» SCA, SKA, SKB , SKC, SSA 


Purpose: These subroutines will compute the heat storage terms for thermal 

model elements of either rectangular, cylindrical, conical or 
spherical geometry. The subroutines are identical except for the 
element voliime . 


Equation: The heat storage term added to the heat balance is: 

« = ifi O S '' ® 

where V is the element volume (in.^). 


Data Cards : Symbolic Name 

NN 
NM 

X, Y, Z, 6) 

R0 , RI , RM 

® j 4’ 1 f 4' p I 

B ^ ^ > 


Value 

node number 

material properties curve relative number 

element dimensions (in. or deg.), 
used to obtain V 
arbitrary constant 


Required Curves: _ 

p - (NM, 1, O) - density (lb /ft^) 

C - (NM, 2, O) - specific heat at constant pressure (Btu/lb -®R)’ 
p m 

Restrictions and Notes: 

1) All subroutines are available in general format niimber 2. 
Subroutines SRA, SCA and SKA are also available in general 
format number 1 . 

2) Curve argument number 2 is available as the temperature of 
node NN. 


Subroutine Geometry: Subroutine Geometry Type 


SRA 

Rectangular 

SCA 

Cylindrical 

SKA 

Conical 

SKB 

Coni cal 

SKC 

Conical 

SSA 

Spherical 
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General Heat Transfer Program 

SRA - Heat Storage - Hect. 

- 

‘^SRA ~ ® 

C ^ ^ ^ B (T- 
p 1728 ^ NN 

- T°j^)/(At/3600) 

NN - Node Number Required curves- 

X 

Y^- Element Dimensions p _ (h„^ 0 ) _ Density 

C-(NM, 2 , 0 )- Specific Heat 

- 




B - Constant ^ 

tl 
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General He.t Transfer program SCA - Heat Storage - Cyl. 

n - „ ^ l6_ (R«5 -RI'^)Z 

^SCA p 360 . 1728 

NM - Node Nvunber 
R0 - Outer Radius 
RI - Inner Radius 
3 - Axial Length 

9 - Angle 

B - Constant 



Required Curves: 

P - (NM, 1, 0 ) - Density 

Cp - (NM, 2, 0 ) - Specific Heat 
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General Heal Transfer Ptoiiam SCA - Heat Storage - Cyl. 

O = „ P ”0 (r0^ - RI 

^SCA P3go. 1728 

^)Z 

— ® - T^) /(At/3600) 

NN - Node Number Required Curves; 

M - Outer Radius p - (NM, -1, 0) - Density 

RI - Inner Radius C - {NM, 2, 0) - Specific Heat 


y Z - Axial Length ^ 

0 - Angle 

B - Constant 
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j^EI 


30 

34 

38 

42 

46 

50 

54 

58 

62 

66 

70 

74 

78 

■ 

SCA 


s 



2 














n 



■ 

NCARD 

H 

NN 

18 

22 

NM 

26 

30 

[R0] 
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Geneial Heat Tiansfei Ptoetain SKA - Heat Storage - Con. 

^ 1 lie (Ri2^+ R8l R 

_ Ska -“S 3 360. ^ 

^ (Sm - t°„)/(ax/36oo 

NN - Node Number 
R0 - Outer Radius 
RI - Inner Radius 
Z - Section Height 
0 “ Angle 

B - Constant 

1 1 

) 

Required Curves; 
p - (NM, 1, 0) - Density 

- (NM, 2, 0) - Specific Heat 

mu 
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18 

NPRNT 

22 

NFMT 

26 

30 

34 

38 
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42 

46 
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• 74 
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General Heal Tianslef Pragiam SKA - Heat Storage - Con. 

„ . „ 1 u0 (R0^ 

^KA ^ p 3 360 . 

t R5) RI-2 RI^)2 , . 0 > ^ . 

1728 « <'^NN - TSj,)/(^t/3600) 

NN - Node Number Required Curves: 

R0 - Outer Radius p - (nm, 1, 0) - Density 

RI - Inner Radius C - (nM, 2, 0) - Specific Heat 

Z - Section Height ^ 

6 - Angle 

B - Constant 
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Gcneial Heat Ttanslei Piogiam SKB - Heat Storage - Con. 

Q = o C ^ RM -1- RM^-TRI^)2 fT' - T° )/Ut/T600) 

"^SKB ^ p 3 360. 1728 ^ NN ^ ^ 

NN - Node Number Required Curves: 

R0 - Outer Radius p - (NM, 1, 0) •• Density 

RI - Inner Radius C - (NM, 2, 0) - Specific Heat 

RM - Middle Radius ^ 

^ 2 - Section Height 

^ ® “ Angle 

i ^ ” Constant 
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General Heat Ttansfei Pfogiam SKC - Heat Storage - Con. 



(R0-RI)^/ Z 

* ® ^'^NN ■ TSn)/[^"/3600) 

p - (NM, 1, 0) - Density 

Cp - (NM, 2, 0 ) - Specific Heat 

360 1728 

NN - Node Number 

R0 - Outer Radius 
RI - Inner Radius 
Z - Height of Cone 

0 - Angle 

^ ^ ^ 6 - Skin Thickness 
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General Heal Tianslei Ptoj/am SSA _ Heat Storage - Sph. 

_®SSA 

^ ^0 (R0^-RI^) 

= ° 3 i80 1728 '■^NN - T°jj)/(At/3600) 

NN - Node Humber Req^uired Curves: 

R0 - Outer Radius p - (m, 1, O) - Density 

RI - Inner Radius C - (NM, 2, O) - Specific Heat 

0 - Longitudinal Angle P 

\ \ ~ ^stitudinal Angle 

<^3 ^ NN \ 02 “ Latitudinal Angle 

» \/\\ ’ B — Constant 
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A. 2.0 Heat Conduction Subroutines 


Section Page 


A. 2.1 DRA, DCA, DCB, DCC, DKA, DKB, DKC, DKE, DSA, DSB, DSC A-15 

(two nodes, same material) 

A.2.2 DRB A-31 

(two nodes or node-source/sink, two materials, 
interface conductance) 

A. 2. 3 DRC, DCG, DCH, DCl A-34 

(two nodes, interface conductance) 
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DRA, DCA, DCB, DCC, DKA. DKB, DKC, DKE, PSA, DSB, DSC 


Purpose ; 


Equation : 


These subroutines will compute the heat conduction terms between 
two elements of the same material. The thermal model elements 
may be of either rectangular, cylindrical, conical or spherical 
geometry. The subroutines are identical except for the conduction 
geometrical shape factor. 


The heat conduction term added to the heat balance is : 



k (A/L) B (Tj^p 


■T 

NN 


) 


where (A/L) is the element conduction geometrical shape factor (in.). 


Data Cards : Symbolic Name 

m 

NP 

NM 

X, Y, Z, (S 
R0, RI, RM 

e , <(>2 

D 


Value 

node number of one element 

node number of other element 

thermal conductivity curve relative number 

element dimensions (in. or deg.) 
used to obtain (A/L) 
arbitrary constant 


Required Curves ; 

k - (NM, 3, O) - thermal conductivity (Btu/hr-ft-°R) 

Restrictions and Notes; 

1) All subroutines are available in general format number 2 . 
Subroutines DRA, DCA, DCB and DCC are also available in 
general format number 1. 

2) Curve argument number 2 is available as the average temperature 
of the two nodes NN and NP. 

3) The heat conduction terms are stored as two-way terms in the 
heat balance. 

U) The subroutines are not valid for METH0D = 7. 


Subroutine Geometry: 


1 routine 

Geometry Type 

DRA 

Rectangular 

DCA 

Cylindrical 

DCB 

Cylindrical 

DCC 

Cylindrical 

DKA 

Conical 

DKB 

Conical 

DKC 

Conical 

DKE 

Conical 

DSA 

Spherical 

DSB 

Spherical 

DSC 

Spherical 
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General Heat Transfei Ptogram DRA - Conduction, One Material - Rect. 

. . ... 

_%RA = ^ ^ 


2 X. 

NN - Node Number Required Curves: 

NP - Node Number k - (NM,3,0) - Thermal Conductivity 

X ) 

Y , - Element Dimensions 
Z ) 

B - Constant 

- 


SsS 

B 

u 

liHii 

m 
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NCARO 
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NPRNT 

22 

NFMT 

26 

30 

34 

38 

42 

46 
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62 
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70 

74 
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Bi 

DRA 
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J 
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L 












■ 
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NCARD 
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NN, 

18 

WP 
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26 
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[X] 

34 
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[B] 
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it6 

R0 + RI 







General Heat Jransler Program 


Z-^ 2_ B (T - T ) 

+ RI 12 ' NP NN' 

NN - Node Number 
NP - Node Number 

__ R0 - Outer Radius 

RI - Inner Radius 

i^Ty 2 - Axial Length 

0 - Angle 

B - Constant 


DCA, Heat Conduction - Cyl., Angular 


Rerjuired Curves : 

k - (NM,3,0) - Thermal Conductivity 


RI, R0, e 9< 0 
RI 5* R0 


NPRNT NFMT 

18 22 26 30 34 38 42 46 50 54 58 62 66 7 0 7 4 78 




NN NP NM 
18 22 26 30 


54l 58 62 66 70 74 78 80 
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General Heal Tiansfer Program DCB - Heat Conduction - Cyl., Axial 

- <'*W^' 

■^': - 

i™)/12Z 

NN - Node Number 
NP - Node Number 
R0 - Outer Radius 
RI - Inner Radius 
Z - Axial Length 
8 - Angle 

B - Constant 



Required Curves : 

k - (NM,3,0) - Thermal Conductivity 
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General Heal Tiansfei piOEram DCB ■ Heat Conduction - Cyl. , Axial 

_S3CB ^ 360 ^ 

ihi - i.'ode Niimber Required Curves: 

liP - iiode iiujnber k - (NM,3,0) - thermal conductivity 

— R0 - Outer Radius 

~ Inner Radius 

— ^ ~ Axial Length 

e - Angle 

\ \ / — / 3 - Constant 


RC0OT 
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NCARD 
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NPRNT 

22 
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General Heat Tianslei Program 


See ^ l 80 12 ^ ~ Sn ^ ^ ^ 


<.,V‘’V • 7 


NN - Node Number 
NP - Node Number 
R0 - Outer Radius 
RI - Inner Radius 
Z - Axial Length 
0 - Angle 

B - Constant 


Dec - Heat Conduction - Cyl., Radial 


Required Curves : 

k - (NM,3,0) - Thermal conductivity 


Note : 

RI, R0 # 0 
RI ^ R0 



I8l 22 


30l 34| 381 4 2 46 50l M 58 6 2 6 6 7 0 7 4 7 8 



NN NP NM 

18 22 26 30 



RI 3 e B 

40 41 M51 60 61 70 71 80 
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General Heat Transfer Progiam DCC - Heat Conduction - Cyl., Radial 

_%CC 

^ H fT T 

l80 12 ^ NP ” N1 

)/ln(R(0/Rl) 

HN - Mode Nximber Required Curves: 

NP - Mode Number k - (NM,3,0), thermal conductivity 

R0 - Outer Radius 
RI - Inner Radius 
Z - Axial Length 
S - Angle 
B - Constant 


NiV — 

Rq *I ' 1 Mote : 

Rl, m ^ 0 

RI ¥ R0 
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General Heat Transfer Program DKA - Heat Conduction, Con., Angular 

Q 

DKA 

^ if 

(1 - In If) - R0 (1 - in g)] B (Tjjp - /{R0 - RM) 

NN - Node Number Required Curves: 

NP - Node Number , » 

R«! - Outer Radius ^ ~ ^””,3,0) - Thermal Conductivity 

RI - Inner Radius 
2 - Axial Length 

0 - Angle 

V B - Constant 

- 


rr 


"V ^ ^ RM - Outer Radius ° ‘ 

\R ' ; "l *M (smaller) RI , R0, RM, 0 ?! 0 

RC0OT 
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DKA . 
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_ 0 



2j 
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General Heal Transfer piogiam DKB - Heat Conduction - Con., Axial 

®-DKB ^ ik 

RI ( R0 - Rm) 

® (■'np - 

fode Number Required Curves : 

ode Number 

uter Radius k - (NM,3,0) - Thermal Conductivity 

nner Radius 

[R0 - RI RM + RI1 
Lr0 + RI RM - RI. 

NN - 1 
NP - 1 
R0 - C 
RI - I 

/ ..'5^ 2 ■ Axial Length 

\ ® " Angle Note: 

Y — Wr o ° " ^oribx,aiu, ^ R0 ^ 0 

_ "V * V ** ■ Outer Radius RI ?! RM ?! R0 

(smallest) 

RC0DT 
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■ 

NCARD 

14 

18 

NPRNT 

22 

NFMT 

26 

3Q 

34 
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50 
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62 
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74 
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80 

DKB 


0 



2 
















RC0OT 

1 

ill 


NN 

18 

MP 

22 

NM 

26 

30 

[R01 

34 

li 

(a] 

42 

(0l 

46 

[B] 

50 


58 

62 

66 

70 

74 

78 

80 

DKB 












































































































































































— 











































A-25 


MCOON NEL L OOC/GLAS COftr*C»f9ATIOM 









































REPORT MDC A3791 


FASTEMP 


»■ 


FASTEMP 

Vc = TBo 


General Heat Tiansfer Program 


DKC - Heat Conduction - Con., Radial 


a /12 

In [(R(a + RM)/(2 RI) 


B (T - T ) 
' NP 



NN - Node Number 
NP - Node Number 
- Outer Radius 
RI - Inner Radius 
a - Axial Length 
0 - Angle 

B - Constant 
RM - Outer Radius 
( smallest) 


Required Curves: 

k - (NM,3,0) - Thermal Conductivity 


Note : 

RI, RM, R0 / 0 
RI / RM / R0 
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Geneial Heat Ttanslei Piogiam DKE - Heat Conduction - Con., Axial 

r . . 

^ ^ ^(R^ - RI) , ^ 

ode Number Required Curves: 

ode Number 

arger Outside Radius k - (NM,3,0) - Thermal Conductivity- 

mailer Outside Radius 


•DKE 

S 

cos Y in 

COS y-6\ 

cos y-&l 

NN - N 
NP - N 
R0 - L 
RI - S 

^^Wx=-iNN g - Height 

R ’ Z, RI, R0 7 0 

\ ^ “ Skin Tnickness 

- — - — Y - arctan [(R0 - 
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General Heal Transler Program 


DSA-Heat Conduction-Sph . .Latitude 


NN - Node Number 
NP - Node Number 
R0 - Outer Radius 
RI - Inner Radius 

0 - Longitudinal angle 

01 - Latitudinal angle 

02 - Latitudinal angle 

3 - Constant 


Req,uired Curves: 

k-(NM,3,0) - Thermal conductivity 


<«, -II 


NCARO NPRNT 

14 18 22 


Note : 

0 ^ , 02 # 

e 0 

0 , 4^02 


30 1 34| 38 1 42 1 46 l 50j 54P 5^ 62 1 66l 70l 74 78 80 


NN NP NM 

18 22 26 30 


[e] I [0i]| [02] |[B] 1 

42 46 ^ 58 62 66 70 74 78 80 
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General Heal Tfsnsfer program DSB-Heat Conduction - Sph., Longitude 

_^DSE 

- k B (T^p-T^„)/ln[tan(|^^)/tan(|,J^)] 

NN - Node number Required curves : 

NP - Node number 

R0 - Outer radius k-{NM,3,0) - Thermal conductivity 

• RI - Inner radius 

0 - Longitudinal angle 

” Latitudinal angle 
^ \ 02 ” Latitudinal angle 

^ \ ^ \ B “ Constant 

1 1 1 Note : 

0, , 0„ 0 or 180 

0, ^ 02 
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General Heat Ttansfci Progiam DSC-Heat Conduction-Sph . , Radial 

S)SC 

- k — 

“ 180 

"1 

iffiTO 

NN .- Node Number Required Curves: 

NP - Node Number 

R0 - Outer Radius k-(NM,3,0) - Thermal conductivity 

RI - Inner Radius 
0 - Longitudinal angle 

~ Latitudinal angle 
02 “ Latitudinal angle 
\ ' B - Constant 

j \ v*« - 

1 Note: 

. RI . R0 0 
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A. 2. 2 DRB 


Purpose: This subroutine will compute the heat conduction terms between 

two elements or an element ajid a source/sink. Conduction is 
through two materials in a rectangular geometry with an optional 
interface conductance. Terms may be one-way or two-way in the 
heat balance. 


Equation: The heat conduction term added .to the heat balance is: 

Q = (B)(BB)(Tj^p-T^) 

12 Xi ^ 12 X? ^ lUlt 
kiAi ^ 2^2 CAp2 


Data Cards: Symbolic Name 

m 

NP 


NMl 

NM2 

XI 

A1 

X2 

A2 

A12 

NC 

B, BB 


Value 

node nximber of one element (if connection 
is one-way, flow is from this node) 
node number of other element or source/ 
sink curve relative number 

a) if 0<NP£.U000, NP is a node number 
with a two-way connection. 

b) if 000 <: NPi 5000, (WP-1^000) is the 
the curve relative number of the 
source/sink temperature. 

c) if NP> 5000, (NP-5000) is a node 
number with a one-way connection 
from node NN 

thermal conductivity curve relative 

number for the element NN 

thermal conductivity curve relative 

number for the element NP 

element length of node NN (in.) 

element cross-sectional area of node NN(in^) 

element length of node NP (in.) 

element cross-sectional area of node NP (in'^) 

cross sectional area for interface 

conductance between node NN and NP (in'^) 

interface conductance curve relative 

number 

arbitrary constants 


Required Curves: 

kq - (NMq, 3, 0) - thermal conductivity (Btu/hr-ft-°R) 
k2 - (NM2, 3, 0) - therma] conductivity ( Btu/hr-ft-°R ) 

C - (NC, 5, O) - interf.ace conductance (Btu/hr-ft2-°R) if NC > 0 
■ T^p - ([NP-UOOO], 13, ) - source/sink temperature (°R) if UOOO < NPi 

5000 
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Restrictions and Notes: 

1) This subroutine is available only in general format number 2. 

2) Curve, argument number 2 is available as : 

a) the temperature of node NN for the evaluation of , 

b) the temperature of node NP (or sourc^sink) for the 
evaluation of k2, 

c) the average temperature of the two nodes NN and NP 
(or source/sink) for the evaluation of C. 

3) The terms B and BB may not be zero 

U) The subroutine is not valid for METHOD = 7. 

Subroutine Geometry: Rectangular 
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General Heat Transfer Program 


SrB = '"'nP - X^/lc^A^ . iJrVC 




NN - Node Number 

NP - £ I 4 OOO Node Number 

- >kOOO, Relative No. for 
Curve 

Xi - Path Length for Node NN 
A]^ - Conduction Area for Node NN 
X 2 - Path Length for Node NP 
Ag - Conduction Area for Node NP 
Aj^g- Conductance area (if needed) 
B - Constant 
BB - Constant 


Required Curves: ~ 

kq - (NM1,3,0)- Thermal conductivit; 
for Node NN 

kg - (NM2,3,0)- Thermal conductivit; 

for Node NP 
if NC > 0 

C - (NC,5,0)- Interface conductance 
if NP > UOOO 

Tjjp- (NP-1|000 ,13 ,1* )- Source/sink 

temperature " 

Notes ; 

Al> Aj, A]^2> BB, ki, kg, C ?f 0 
Xi and Xg 0 
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A. 2. 3 PRC, DCG, DCH, PCI 


Purpose: These subroutines will compute the heat conduction terms for 

interface conductance between two elements. The thermal model 
elements may be of either rectangular or cylindrical geometry. 
The subroutines are identical except for the interface area. 

Equation: The heat conduction term added to the heat balance is: 

where A is the interface area (in ). 


Data Cards: 


Symbolic Name 


NN 

NP 

NC 

Y, Z 
R2), RI, 
B 


I 

9( 


Value 

node number of one element 

node number of other element 

interface conductance curve relative number 

element dimensions (in. or deg.) 

to obtain A 

arbitrary constant 


Required Curves : 

C - (NC, 5,0)- interface conductance (Btu/hr-ft^-°R) 

Restrictions and Notes: 

1) All subroutines are available in general format numbers 1 and 2.. 

2) Curve argument ntimber 2 is available as the average temperature- 
of the two nodes NN and NP. 

3) The heat conduction terms are stored as two-way terms in the 
heat balance 

These subroutines ' are not valid for METH0D T. 


R r c '.:t i r.e Oe omet ry : Subroutine 

PRC 

DCG 

DCH 

PCI 


Geometry 
Rectangular 
. Cylindrical 
Cylindrical 
Cylindrical 
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General Heat Transfet Ptog(an) DRC - Interface Conductance - Rect. 

Y 2 

3>rc ■ ^ Tinr 


NN - Node Number Required Curves: 

NP - Node Number C - (NC, 5, O)- Interface Conductance 

- Element Dimensions 
■B - Constant 

- 

NN< iNP 



A 
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General Heat Transfer Progiam DRC - Interface Conductance - Rect. 

5 

-DRC ^ litit ® *'^NP ■ 

NN “ Node Number Required Curves: 

NP - Node Number C - (NC,5,0) - Interface Conductance 

y) 

gj Element Dimensions 

^ ” Constant 


NP 
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General Heal Ttansfei Piogram DCG - Interface Conductance-Cyl. , Angular 

QdCg = C B (T,p - 

NN - Node Number Required Curves : 

NP - Node Number 

- R0 - Outer Radius C - (NC, 5,0)- Interface Conductance 

BI - Inner Radius 
® Length 

\ ^ L “ Constant 


L __ _ _ 
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NPRNT 
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NFMT 

26 

30 

34 

38 

42 

46 

50 

54 

58 

62 

66 

70 

74 

78 

80 

Idcg 


0 
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General Heat Transfer Ptogram DCG - Interface Conductance-Cyl. , Angular 

_ *^DCG ” 

„ (R(i5-R 
^ l!. 

T — B (T - T ) 

4 ^ NP NN'^ 

NN - Node Number Required Curves: 

NP - Node Humber 

s=^ R0 - Outer Radius C - (NC, 5,0)- Interface Conductance 

— 7 RI - Inner Radius 

/ / Z - Axial Length 

/ . B - Constant 

Nish »NP 

''' 
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Ijkg __ 


m 



2 






J 









n 


RC(30T 

1 

■ 


NN 

18 

NP 

22 

NC 

26 

30 


[RI] 

38 

[2] 

42 

46 

[B] 

50 

54 

58 

62 

66 

70 

74 

78 

80 

DCG 

























































m 





















m 







































































































































A- 38 


MCOOMtSI£LL OOUGLMS COt9§»OftATiOM 









































REPORT MDC A3791 


FASTEMP 


FASTEMP 

r ~ — — , , , 

Genetal Heat Tiansfei pfog/am DCH - Interface Conductance - Cyl , Axial 

Q = 

■DCH 

„ RI^) 

sis' ® ^'^NP ■ ’’’nn^ 

NN - Node Number Required Curves; 

NP - Node Number 

R0 - Outer Radius C-fNC.5.0) - Interface Conductance 


in; 



1 

0- \ 

. / \ 

o 

NNirNP J RI - Inner Radius 

^ 0 - Angle 

Vs-?' 
V, ; 



B - Constant 
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General Heat Transfer Program 

DCH - Interface Conductance - Cyl, Axial 

“^DCH 

""•A" 
. ^ 

_ „ ti(R0^ - RI^) 0 , V 

■ ^ nn; — w ® 

NN - Node Number Required Curves: 

NP - Node Number 

7 R0 - Outer Radius C-(NC,5,0) - Interface Conductance 

! NNVNP^ RI - Inner Radius 

\ 0 - Angle 

' B - Constant 

r* o ^ 
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General Heat Translei Pfogtam DCI - Interface Conductance - Cyl. , Radial 

- _2tiR0Z 0 pfrT, 

%CJ IUI 4 3 S 0 ® 

NN - Node Number Required Curves: 

NP - Node Number C-(NC,5,0) - Interface Conductance) 

'V ^ 

%> \-eJ * 

- jS--- 

Kfi - uuter naaius 
Z - Axial Length 
0 - Angle 

B - Constant 
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General Heat Ttansfei Piogiam DCI _ Interface Conductance - Cyl. , Radial 

2ti R0 a 0 , , 

^DCI ^ Ihh 360 

NN - Node Number Required Curves : 

NP - Node Number 

j — / Z - Axial Length 

\ Q - Angle 

^ B - Constant 

RC0OT 
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A. 3.0 Heat Convection Subroutines 


Section 


Page 


A. 3.1 VRA, VGA, VKA, VSA A-44 

(forced convection heat transfer coefficient, 
adiabatic wall temperature) 
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A. 3.1 VRA. VGA, VKA, VS A 


Purpose: These subroutines will compute the heat convection terms for 

thermal model elements of either rectangular, cylindrical, conical 
or spherical geometry from a source/sink. The subroutines are 
identical except for the element surface area. Heat transfer coeffi- 
cient is provided as known function. 

Equation: The heat convection term added to the heat balance is: 

Q = 1 h A B (T -T„„) 

^ 144 ' aw NN 

where A is the element surface area (in^) 


Data Cards : Symbolic Name 

M 
NH 

NTAW 

Y., Z 
Re), Ri 

<t>2 

I 

B 


Value 

node number 

heat transfer coefficient curve relative 
number 

adiabatic wall temperature curve relative 
nvur'ber 

element dimensions (in. or deg.) 

area option (VCA only) 
arbitrary constant 


Required Curves : 

h - (NH, 15, O) - heat transfer coefficient (Btu/hr-ft^-°R) 

T - (NTAW, 13, l) - adiabatic wall temperature (°R) 
aw 

Restrictions and Notes: 

1) All subroutines are available in general format number 2. 
Subroutines VRA and VKA are also available in general format 
number 1 . 

2) Curve argument number 2 is available as the temperature of 
node NN. 

3) These subroutines can be used to force a node to follow a 
prescribed temperature by letting the T^^ curve be the 
prescribed value and making the heat transfer coefficient 
a large number (e.g. 10°). 


Subroutine Ceometry': Subroutine 

VRA 

VCA 

VKA 

VSA 


Geometry 

Rectangular 

Cylindrical 

Conical 

Spl’erical 
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Geneial Heal Transfef Proc(am " Convection Using h and Taw Curves - 

Rect. 

5 

Y ? 

= h -Vr B (T - T 
VRA II 4 I* AW 

/ NN 
Y / • 

// / 

NN ^ 

NN - Node Number Required Curves ; 

“ Element Dimensions ^ “ (NH,15,0) - Heat transfer 

* — / B - Constant coefficient 

^ - (NTAW,13,1) - Adiabatic wall 

/' ] temperature 
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General HeaiTranslei Program ™ " Convection Using h and Taw Curves - 

Rect. 

®VRA Ihk ® 

NN — Node Number Required Curves : 

Y) 

2 > - Element Dimensions . ,c n\ u ^ a. 

, A h - (NH,15,0; - Heat transfer 


/ / 

B - Constant coefficient 

*^AW ~ (NTAW,13,1) - Adiabatic wall 
temperature 

- 

Y 

NN 

/ • 

/ 



/ 

/ 



■ 
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General Heal T.ansle, Program " Corrvection Usirrg h and Taw Curves - 

Cyl. 

Q. 

VGA ■ unr ® ^'^AW”'^NN^ 

^1-0 

NH - Node Number Required Curves: 

R0 - Outer Radius . c„u ^ ^ 

h - (NH,15,0) - Heat transfer 
RI - Inner Radius ^ 

„ , . I T coefficient 

e ~ An^le '^AW ~ (»TAW,13,1) - Adiabatic wall 

B IcoLtant temperature 

I - 0, use Area = -jgQ 

1, use Area = (R0^-RI^) 

j3/^^ 

RCjSOT 
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G6nerat He 3 t Tr 3 nsfer pfogram VKA - Convection Using h and Taw Curves 

- Con, 

Q - J, ^ iRg-Ri) Va^ > (R{)-R] 

JVKA W lUU 

NN 

< V — 

NH - Node Nmnber 

R(8 - Outer Radius 

RI - Inner Radius 

Z - Conic Section Height 

C - Angle 

B - Constant 



Required Curves: 
h - (NH,15,0) - Heat transfer 
coefficient 

'^AW ~ (NTAW,13,1) - Adiabatic wall 
temperature 
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0 = h LR0-. B lJ.Ygi B (T T ) 

NN - Node Number 
R0 - Outer Radius 
RI - Inner Radius 


'vv / ^ 



Geneial Heat Transter Progiam 


VKA - Convection Using h and Tav Curves - Con. 


Required Curves: 
h - (NH,15,0) - Heat transfer 
coefficient 


Z - Conic Section Height - (NTAW,13,l) - Adiabatic wall 


e - Angle 
B - Constant 


temperature 


NPRNT NFMT 
18 22 26 



30 3 4 38 42 46 50 5 4 58 6 2 66 7 0 7 4 78 




181 22 


[R0] [RI] 
30 34 38 


[e] [B] J 

46 50 54I 58 62 66 70 74 78 
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General Heat Transler Program 


VbA - i-onvection Using h and Taw Curves 

- S 


3vsa = 180 IW ® 


' -T ) 
AW NH' 



NN - Node Number Required Curves : 

R0 - Outer Radius h - (NH,15,0) - Heat transfer 

8 - Longitudinal Angle coefficient 

01 - Latitudinal Angle Taw - (NTAW,13,l) - Adiabatic 

02 - Latitudinal Angle wall temperature 

B - Constant 
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A. 4.0 Heat Radiation Subroutines 


Section Page 


A. 4.1 

RRA, RCA, RKA, RSA 
(external radiation, eF) 

A-52 

A. 4. 2 

RRB 

(internal radiation, two gray plates) 

A-59 

A. 4. 3 

RRC, RCC, RKC, RSC 
(internal radiation, eF) 

A-62 

A. 4. 4 

RRD 

(internal radiation, enclosure) 

A-69 
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A . 4 . 1 RRA, R C A. RFCA, RSA 


Purpose: These subroutines will compute the heat radiation terms to thermal 

model elements of either rectangular , cylindrical, conical or 
spherical geometry from a source/sink. The subroutines are 
identical except for the element surface area. 


Equation: The heat radiation term added to the heat balance is: 

Q = iw ^ F A ® (V-Tp,j^M 

Q 

where a is the Stefan-Boltzmann constant (0.171^ x 10” Btu/hr- 
ft2_°R^) and A is the element surface area (in^). 


Data Cards : Symbolic Name 

NN 
NF 


NM 

NTR 


X , X. 

R0, RI 

4'i 5 <t>2 
I 
B 


Value 

node number 

shape factor curve relative number from 
node M to source/sink 
eipissivity curve relative number 
radiation source/sink temperature 
curve relative number 

element dimensions (in. or deg.) 
used to obtain A 
area option (RCA only) 
arbitrary constant 


Required Curves : 

F - (NF, IT. 0) - shape factor (-) 
e - (NM, 4, 0) - emissivity (-) 

Tr - (NTR, 13, 2) - radiation source/sink temperature (°R) 


Restrictions and Notes: 

1) All subroutines are available in general format number 2. 
Subroutines RRA and RKA are also available in general 
format number 1. 

2) C\irve arg\mient number 2 is available as the temperature of node NN. 

3) These subroutines are not valid for METH0D = 5- 

4) Approximations are made for T^ in order to retain a linear 

set of equations. These approximations are given in Appendix B. 
These approximations require that the temperature change 
of the radiation nodes remain small over the transient 
time steps. For steady state solutions the approximation is 
exact . 


Subroutine Geometry: 


) routine 

Geometry 

RRA 

Rectangular 

RCA 

Cylindrical 

RKA 

Conical 

RSA 

Spherical 
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General Heat Ttanslei Piogtam RRA - External Radiation - Rect. 

Q 

RRA = ° ^ ^ M ® 

NN 

Y / • 

/ 

- t' ) 

NN - Node number Required curves: 

Yl F - (NF,17,0) - view factor 

gj- -Element Di^nsions c - (NM,4,0) - emissivity 

^ B - Constant Tr - (NTR,13,2) source/sink 

* / temperature 

1 / 
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General Heat Transfer Program 


RCA - External Radiation - Cyl, 





/ < \- 
V. V \ 


\ / 


NN - Node number 
R0 - Outer radius 
RI - Inner radius 
Z - Axial length 
0 - Angle 

B - Constant 


Required curves ; 

F - (HF,17,0) - view factor 
E - (NM,ir,0) - emissivity 

T^- (NTR,13,2) - source/sink 
temperature 


I - 0, use Area - R^Z Note: Approximation used for T 

T A rr9 / ^2\ 

1, use Area = (R^-R^) 



NPRNT NFMT 
18 22 26 
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General Heal Tfansfet Pfogtam rSA - External Radiation - Sph.. 

Asj 

i = ° " ilo fw ® 

NN “ Node number Required curves: 

R0 - Outer radius F - (NF.IT.O) - View factor 

0 - Longitudinal angle e - (NM,U,0) - Qnissivity 

' 0^ - Latitudinal angle T - (NTR,13,2) - Source/sink temp. 

\ 0? ■ Latitudinal angle 

\ ® ” Constant 

\ % U Mote; Approximation used for T** 
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A. 4. 2 


Purpose : 


Equation; 


Data Cards : 


Required Curves : 

e, - (NMl, 4, 0) - emissivity {-) 

E 2 “ (NM2, 4, O) - emissivity (-) 

Restrictions and Notes: 

1) This subroutine is available in general format numbers 1 and 2. 

2) Curve argument number 2 is available as: 

a) the temperature of node NN for the evaluation of 

b) the temperature of node KP for the evaluation of Z 2 > 

3) This subroutine is not valid for METH0D = 5, T- 

4) Approximations are made for in order to retain a linear 
set of equations. These approximations are given in 
Appendix B. These approximations require that the temperature 
change of the radiation nodes remain small over the transient 
time steps. For steady state solutions the approximations are 
exact. 

5 ) The heat radiation terms are stored as two-vay terms in the 
heat balance. 

Subroutine Geometry: Rectangular 


Tills subroutine will compute the heat radiation terms between two 
thermal model elements of rectangular geometry. The radiation 
shape factor is computed for two gray infinite parallel plates. 

The heat radiation term added to the heat balance is: 


Q = 


TTHT o Y Z B 


(T, 


NP 


NN 


- 1 ) 

^1 ^2 


where o is the 
(0.1714 X 10“' 


Stefan-Boltzmann constant 
Btu/hr-ft2 - Ob'*). 


Symbolic Name 

NN 

NP 

NMl 

NM2 

Y. a 

B 


Value 

node number of one node 

node number of other node 

emissivity curve relative number of node NN 

emissivity curve relative number of node NP 

element dimensions (in.) 

arbitrary constant 
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General Heat Transfer program pj ^3 _ Internal Radiation - Beet. 

®RRB ° ° 

(-1 * V 1 

cS 

^ B (t‘‘ T** ) 

iW 

NN — Node number Required curves: 

NP - Node number - (NM1,U,0) emissivity for node 

N - Element Dimensions surface- NH 

) “ (^2,4,0) - emissivity for node 

B - Constant ^ surface NP 

h 

— Z Note: Approximation used for T 

Y 








_ 

10 


18 

NPRNT 

22 

NFMT 

26 

30 

34 

38 

42 

46 

50 

■ 

58 

62 

66 

70 

74 

78 

80 

IsEB 


0 



1 







■1 









RC0OT 

1 

■ 

NCARD 

NN 

18 

NP 

22 

NMl 

26 

W!2 

30 

31 40 

Hi 

E 

51 60 

B 

61 70 

71 80 

RRB 


































































































































A-60 


MCDONNmi. L DoucLMm conr*ortA noN 
































FASTEMP 


REPORT MDC A3791 



%RB ° I 1 ^ 1 _ ^ 


General Heal Transfer Program 


M “ <^p - 




Nil - Node mamber 
NP - Node number 

Y) 

g? - Element Dimensions 
B - Constant 


RRB - Internal Radiation - Rect. 


Required curves: 

- {NM1,1*,0) - emissivity for node 
surface NN 

- (NM2,1*,0) - emissivity for node 
surface NP 


Note: Approximation used for 
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Purpose: These subroutines will cOTipute the heat radiation terms between 

thermal model elements of either rectangular, cylindrical, conical 
or spherical geometry. The subroutines 6u*e identical except for 
the element surface area. 

Equation: The heat radiation term added to the heat balance is: 

Q = » c P A B (T^'' . 

Q 

where a is the Stefan-Boltzmann constant ( 0.171^ x 10”° Btu/hr - 
ft2 _ or 4) and A is the element surface area (in^). 

Value 

node niimber of one element 
node number of other element 
emissivity curve relative number 
shape factor curve relative number from 
node NN to node NP 

element dimensions of node NK (in. or deg.) 

area option (RCC only) 
arbitrary constant 

Required Curves: 

e - (NM, U, 0) - emissivity (-) 

F - (NFA, 17, 0) - shape factor (-) 

Restrictions and Notes: 

1) All subroutines are available in general format number 2. 
Subroutines RRC and RKC are also available in general format 
number 1. 

2) Curve argument number 2 is available as: 

a) the temperature of node NN for the evaluation of e, 

b) the temperature of node NP for the evaluation of F. 

3) These subroutines are not valid for METH0D =5,7 

1^) Approximations are made for T^ in order to retain a linear set 
of equations. These approximations are given in Appendix B. 

These approximations require that the temperature change of 
the radiation nodes remain small over the transient time steps. 
For steady state solutions the approximations are exact. 

5) The heat radiation terms are stored as two-way terms in the 
heat balance. 

Subroutine Geometry: Subroutine 

RRC 
RCC 
RKC 
RSC 


Geometry 

Rectangular 

Cylindrical 

Conical 

Spherical 




A-62 

MCOOMMCLC. DOUGLAa COmfOfW AXIOM 



FASTEMP 


REPORT MDC A3791 


FASTEMP 

Genctal Heat Tianslei Ptogram RRC - Internal Radiation - Rectangular 

c 

(T** - T** ) 

NP NH' 

NN - Node number 
NP - Node number 

H - 

2f Element Dimensions 

B - Constant 

1_ 

Required curves; 

F - (NFA,17iO) - Vie« factor 
E - - Quissivity 

Note; Approximation used for T 
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FASTEMP 

General Heal Transfer Program RRC - Internal Radiation - Rectangular 

3rrc " ° 

^ ieI “ <■'» 

— z- 

NN - Node number Required curves: 

NP - Node number F - (NFA,17,0) - view factor 

. £ - (NM,i,0) - emissivity 

21 ~ Element Dimensions 

B - Constant 

Note: Approximation used for T** 
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General Heat Tfanster program 


RCC - Internal Radiation-Cyl. 




I-JN - Node number 

NP - Node number 

R0 - Outer radius node NN 

RI - Inner radius node NN 

2 - Axial length 

6 - Angie 


b - Constant 

I - 0, use area = ^ R^Z 



Required curves: 

F - (UFA, 17,0) - View factor 
c - (NM,U,0) - Emissivity 

Note: Approximation used for 
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Gencol Heal T^nsfe; Ptogiam RKC - Internal Radiation - Con. 

T,e (R0 ^ Ri)y[i 
^KC 360 ] 

0 ) 

^ + (R0 - P.I)^' ,„k U 

Ui. ■ ^ 

NN - Node niunber 

NP - Node number 

R0 - Larger radius node NN 

RI - Smaller radius node NN 

Z - Conical section height 

0 - Angle 

B - Constant 



Required curves: 

F - (NFA,1T,0) - View factor 
e - (NM,L,0) - Emissivity 

Note: Approximation used for 
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FASTEMP 

General Hea! Ttanslei Pragiaiii 

RKC - Internal Radiation - Con. 

„ u8 

^KC ■ ° ^ 360 

(R0 + RI) Vz^ + (R0 - RI)^ 

(^!p - 

e number Required curves: 

e number F - (NFA, 17,0) - view factor 

ger radius node HN e - (KM,1*,0) - emissivity 

Her radius node NN 

ical section height Note: Approximation used for 
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General Heat Transtet Ptogram RSC - Internal Radiation - Sph. 

Ssc = ° ^ llo ife ® 

NN - Node nujnber Req^uired curves : 

NP - Node number F - (NFA,17,0) - View factor 

/ ^ R0 - Outer radius node NW c - (NM,i,0) - Qnissivity 

- Smallest latitudinal 

V. angle Note: Approximation used for 

^ \ " Lsirgest latitudinal 

_ \ \ ^ B - Constant _ 

y >-■ 0 
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A. 4. 4 RRD 


Purpose: This subroutine will compute the heat radiation terms between a 

thermal model element and up to three other elements in an 
enclosure. 


Equation: The heat radiation terms added to the heat balance are: 

where c is the Stefan-Boltzmann constant 
(O.lTlU X 10-8 Btu/hr - ft^ - °R^). 


Data Cards: Symbolic Name 

NN 

Nl, N2, N3 
A 

FNl, FN2, FN3 


Required Curves : 
None 


Value 

node number 

node numbers of other surfaces (up to three) 

surface area of node NN (in^) 

shape factors from node NN to node Nl , 

N2 and N3 respectively 


Restrictions and Notes: 

1) This subroutine is available in general format numbers 1 and 2. 

2) This subroutine is not valid for METH0D =5,7. 

3) Approximations are made for T** in order to retain a linear 

set of equations. These approximations are given in Appendix B. 
These approximations require that the temperature change of 
the radiation nodes remain small over the transient time steps. 
For steady state solutions the approximations are exact. 

4) The heat radiation terms are stored as two-way terms in the 
heat balance. 


Subroutine Geometry: Rectangular 
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Gcneial Heal Ttanslei piogtam 


RRD - Internal Radiation - Rect. 


.1* 'U . 

Srd " TITir ■ '^NN 



NN - node number 
N1 - Node number 
N3 - Node number 
N3 - Node number 
*NH“ Area of Node NN 
FNl- View factor node 
NN to N1 

FIJ2- View factor node 
NN to N2 

FN3- View factor node 
NN to N3 


Required curves; 
None 


Note: Approximation used for T 
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General Heal Transler Program 


RRD - Internal Radiation - Rect. 


%RD = ° ^NI ^ 



NN - 

Node 

munber 

N1 - 

Node 

number 

N2 - 

Node 

number 

N3 - 

Node 

Number 

ANN- 

Area 

of Node NN 

FNl- 

View 

factor from 


node 

NN to N1 

FN2- 

View 

factor from 


node 

NN to N2 

FN3- 

View 

factor from 


node 

NN to N3 


Required curves: 
None 


Note; Approximation used for T 
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A. 5.0 Heat Flux Subroutines 


Page 

QCRA, QCCA, QCKA, QCSA A-73 

(heat flux, one multiplier) 
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A. 5.1 QCRA, QCCA, QCKA, QCSA 


Purpose: These subroutines will compute the heat flux terms for thermal 

model elements of either rectangular, cylindrical, conical or 
spherical geometry with a prescribed heat flux imposed on its 
surface. The subroutines are identical except for the element 
surface area. 


Equation : 


The heat flux term added to the heat balance is: 


Q = 


Ihh 


q A P B 


where A is the element surface area (in ), 


Dat a Car ds : Symbolic Name 

NN 

NPER 
NQ 
Y, Z 
R5i, RI 

®>4>i,<(>2 

I 

B 


Value 

node number 

multiplier curve relative number 
heat flux curve relative number 

element dimensions (in. or deg.) 

are i option (QCCA only) 
arbitrary constant 


Required Curves : 

P - (NPER, 99 , O) - multiplier (-) 
q - (NQ, 99 , 0) - heat flux (Btu/hr - Ft^ ) 


Restrictions and Notes: 

1) All subroutines are available in general format number 2. 
Subroutines QCRA and QCKA are also available in general 
format number 1. 

2) Curve argument number 2 is available as the temperature of 
node NN. 


Subroutine Geometry: Sub routine Geometry 


QCRA 

Rectangular 

QCCA 

Cylindrical 

QCKA 

Conical 

QCSA 

Spherical 


ORlOjS'iAL PAGE IS 
OF POOR QUAUTY 
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General Heat Tianslet Piogiam QCRA - Q„,,o„„ - Rect . 

LUKVb 

_^QCRA ^ 

is! ^ ® 

NN - Node number Required curves: 

^ . P - (NPER,99,0) - Multiplier 

y, g/ - --^enent Dimensions q _ (fIQ,99,0) - Heat flux 

/ B - Constant 
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General Heat Transfer program QCCA - Q Curve - Cyl . 

VCA = ^ ^ ^ » 

NN “ Node number Required curves: 

R0 - Outer Radius P - (nPER,99,0) - Multiplier 

^I«0 - Inner Radius q - (NQ,99,0) - Heat flux 

\ 8 - Axial length 

7"1 ^ Q - Angle 


V ^ ~ ^ Area = R0 2 

- 

* X 

V® 1 use Area = ^ (R0^ - RI^) 

<l V / 360 
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General Heat Tiansfet P/osiam qCKA - Q Curve - Con. 

-6 (R0+RI) + {R0-RI 

5iCKA 360 iTT 

NN 

< ^ 

^ PB 

KU - Mode number 
R0 - Larger outside radius 
RI - Smaller outside radius 
E - Conic section height 
0 - Angle 

B - Constant 



Required curves: 

P - (NPER,99,0) - Multiplier 
q - (NQ,99.0) - Heat flux 

L_.. 

■H 

10 


18 

NPRNT 

22 

NFMT 

26 

30 

34 

33 

42 

46 

50 

54 

58 

62 

66 

70 

74 

78 

1 

5CKA 


0 



i 















n 


RC0DT 

1 

■ 


NN 

13 

HP2R 

22 

■ 

30 

Hi 


Hni 

[e] 

61 70 

■BhI 

jcka 

■ 

'm 










Zj 













__| 
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General Heat Transfer progtam QCKA - Q CURVE - Con. 

n6 

(R0 + RI) Va^ + (R0 - RI)^' 

PB 

e number Required curves: 

ger outside radius P - (NPER,99,0) - Multiplier 

Her outside radius q - (NQ,99.0) - Heat flux 


ICKA 360 

NN 



lUU 

NN - Nod 
R0 - Lar 
RI - sma 

-'"V ® - Angle 

^ B - Constant 


1 

NCARO 

14 

18 

NPRNT 

22 

NFMT 

26 

30 

34 

38 

42 

46 

50 

54 

58 

62 

66 

70 

74 

78 

■ 

QCKA 


0 



2 














n 


RC0OT 

1 

■ 

NCARO 

14 

NN 

18 

■ 

NQ 

26 

30 

[R0] 

34 

[RI] 

38 

[2] 

42 

[8] 

46 

[B] 

50 

54 

58 

62 

66 

70 

74 

78 

■ 

QCKA 
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General Heat Tianslet Pfogiam QCSA - Q Curve - Sph. 

^QCSJ 

^ ilo fw 

NW — Node number Required curves : 

R0 - Outer radius P - (NPER,99,0) - Multiplier 

e - Incremental longitudinal q - (NQ,99,0) - Heat flux 
\ ^ angle 

- Smallest latitude angle 
'^o ■ Largest latitude angle 
^ ” Constant 

/> JT ' 




I 

NCARD 

14 

18 

NPRNT 

22 


30 

34 

38 

42 

46 

50 

54 

58 

62 

66 

70 

74 

78 

■ 

QCSA 


0 



2 














■ 



■ 

mu 

NN 

18 



30 

[R0] 

34 

33 

[0] 

42 

[0,1 

46 

[0p] 

50 

[B] 

54 

58 

62 

66 

70 

74 

78 

80 

QCSA 
































1 

1 
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■ 
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A. 6.0 Fluid Flow Subroutine 


Section Page 

A. 6.1 FLUIDS A-81 

(fluid convection with wall heat transfer) 
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Purpose : 


Equations : 


A. 6.1 FLUIDS 

This subroutine will compute the heat terms for a fluid flowing 
by one or more surfaces. The fluid to wall heat transfer coeffi- 
cient h may be input by the user as curve data or may be computed 
by the subroutine. Options are provided in the subroutine for: 
Option 0 - bulk temperature method 
Option 1 - outlet temperature method 
Option 2 - modified h method 
Option 3 - modified m and h method 
This subroutine can also be used to compute the pressure drop 
in the fluid control volume. 

The heat terms added to the heat balance are: 


3600 - T ”) - Q +0 

1728 At F F '' ^m 


where : 


Option 0 - 0^ * m (T^^ - T^^) 

hA B 

%, ^ ' 144 ~ 


T, = 1/2 (T^^ -H T^^) 


Option 1 - Q = m C (T* - T' ) 

m p 3 NI N0 

hA B„ 


0 = (T' - T' ) 

144 ^ NW 


T„ = T 


N0 


Option 2-3 = 


hA B 

(E ^)/A Cp B3 


Q_ - 


m 




hA B 


^F ^N0 
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Option 3 - & = a -- - -) /m C B. 

144 p J 


c = g - - - " 


B(e^ - 1) 


Q^ = ACpB3 ta - C)/(e'^ - DlHT^j - V 


Q„ = 1:^ - C T-^ - (1 - C) t; ] 

144 


~ ^N0 


The fluid element control volume equation options are: 
user specified 

V = XI 
rectangular 

V = XI • X2 * X3 
cylindrical 

V = 7 . X2^ . XI 

4 


If the fluid to wall heat transfer coefficient is computed by the 

subroutine, the equivalent diameter, D j'and the~cross sectional 

e 

area, A^, options are: 

user specified 

D = X2 
e 

A = X3 
c 

rectangular 

= 2 • X2 ’ X3 
e X2 + X3 

A = X2 •X3 
c 

cylindrical 

D = X2 

A = 7 • X2^ 
c 4 

The wall area equation options are: 
user specified 
A = X4 
rectangular 
A = X4 • X5 
cylindrical 
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A = TT • X4 • X5 

If the curve relative number IH is not zero, the heat transfer 
coefficient must be input by the user as curve data. If IH is 
zero, the heat transfer coefficient will be computed by the sub- 
routine according to the relations 


N- = 


EXLl 


Pr 


EXL2 


D 

o 


EXL3 


w 


EXL4 


for the laminar case when Re £ R®/., > 

OIj 

EXT3 EXT4 

»„ = . e, • (^) 

o w 


for the turbulent case when Re > Re 

Nu = logarithmic Interpolation between the laminar Nusselt 

number evaluated at Re and the turbulent Nusselt num- 

wLi 

ber evaluated at Re^^. 

X 

when Re < Re < Re„„, where; 

1 


h D 


Nu = Nusselt Number = 


Re = Reynolds number = 


12 

p, V^D 

b e 


Sb, 


12 m D 


3600 


Pr = Prandtl Number = 


3600 Cp, p, 

O D 




A 

c 


b and w denote properties evaluated at the bulk fluid and wall 
temperatures, respectively. 


If the user specifies a value for ILC in the subroutine zero card, 

the constant and exponents for the laminar Nusselt number equation 

are obtained from the dimension table as follows; 

C, ... = value in dimension table location ILC 
LAM 

EXLl = value in dimension table location ILC+1 

EXL2 = value in dimension table location ILC+2 

EXL3 = value in dimension table location ILC+3 

EXL4 = value in dimension table location ILC+4 


If ILC is zero or is not specified, the following default values 
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will be used in the laminar equation: 

"lam = 

EXLl = 0.3333333 
EXL2 = 0.3333333 
EXL3 = 0.3333333 
EXL4 = 0.140 

If the user specifies a value for ITC in the subroutine zero 
card, the constant and exponents for the turbulent Nusselt number 
equation are obtained from the dimension table as follows: 

Ctittoo “ value in dimension table location ITC 

EXTl = value in dimension table location ITC+1 

EXT2 = value in dimension table location lTC+2 

EXT3 = value in dimension table location lTC+3 

EXT4 = value in dimension table location ITC+4 


If ITC is zero or is not specified, the following default values 
will be used in the turbulent equation: 


^TURB 

EXTl 

EXT 2 

EXT3 


0.0225 

0.80 

0.3333333 

0.0 


EXT4 = 0.140 


If the user specifies a value for IF in the subroutine zero card, 
the pressure drop in the fluid control volume will be computed. 
The equation is: 


AP 


1 . 4fL (m/3600) ^ 

e 2 p(A^/144) 


psi 


P 

out 


where: 



AP 


f = friction factor input as a curve versus Reynolds number 
(curve argument 30) 

L = XI variable 

The inlet pressure is optional and is obtained from the card 1 
type. If the inlet pressure is specified, the outlet pressure 
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replaces it allowing the pressure through a series of control volumes 
to be computed. The inlet pressure to the first control volume can 
be specified by the M0PA subroutine. 

If the pressure drop option is specified, the following printout 
is obtained for each fluid control volume at the print times 
specified by curve 1, 10, 2. 

FLUIDB NI N0 AP P Re Nu f 

out 

(The Nusselt number is without the viscosity factor, 
which is dependent on each wall section.) 


Data Cards ; Symbolic Name 
NCF0L 
IFM 


RECL 

RECT 

ILC 

ITC 

IF 

Card 1 LC 


NI 

N0 

IFP 

IMD 

IGV 


Value 

number of cards to follow 
fluid method 

= 0 bulk temperature method 
= 1 outlet temperature method 
= 2 modified h method 
= 3 modified m and h method 
critical Reynolds number for laminar flow 
critical Reynolds number for turbulent flow 
index for laminar Nusselt number constants 
or zero 

index for turbulent Nusselt number constants 
or zero 

friction factor curve relative number or 
zero 

last card in control volume 
= 0 no 
= 1 yes 

fluid inlet node number 
fluid outlet node number 
fluid property curve relative number 
mass flow rate curve relative number 
control volume geometry option 
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Card 2 


B1 

XI 

X2 

X3 

XL0 

B3 

PIN 

LC 


NW 

IH 

IGA 


B2-- - 

XA I 
X5 f 

Required Curves: 
f 

P 

R 

C 

P 

k 

P 

m 

h 

Restrictions and 
1 ) 

2 ) 


= -1 user specified 
= 0 rectangular 

= 1 cylindrical 

volume term multiplier 

element dimensions for volume 

flow development length 
flow term multiplier 

control volume inlet pressure if IF 0 
last card in control volume 
= 0 no 
= 1 yes 

wall node number 

wall heat transfer coefficient curve 

relative number or zero 

area geometry option 

= -1 user specified 
= 0 rectangular 

= 1 cylindrical 

area term multiplier 
element dimensions for area 

- (IF, 27, 0) - friction factor if IF # 0 

- dlFP|, 1, 0) - fluid density (Ib/ft^) if IFP > 0 

- (|IFP|, 1, 1) - gas constant (ft-lb^/lb^ - °R) if IFP < 0 

- (|IFP|, 2, 0) - fluid specific heat (Btu/lb^ - °R) 

- ( I IFP I , 3, 0) - fluid thermal conductivity (Btu/hr-f t-°R) 

if IH = 0 

- (IiFPI, 8, 0) - fluid viscosity (lb /ft-sec) if IH = 0 

‘ ‘ m 

- (IMD, 19, 1) - fluid mass flow rate (Ib^/hr) 

- (IH, 15, 2) - wall heat transfer coefficient (Btu/hr-f t^-°R) 
if IH 0 

Notes: 

This subroutine is only available in general format 2. 

Curve argument number 2 is available as the bulk fluid tem- 
perature Tg = 1/2 (Tj^j + For evaluation of p^, 
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curve argument number 2 is available as the wall temperature 

3) If the presure drop option is used, curve argument number 30 
is available as the Reynolds number for the fluid control 
volume. 

4) For IFF < 0, fluid density is computed using the equation 

p = P/RT . Pressure is obtained from curve argument number 

D 

12 . 


5) 

6 ) 


7) 


8 ) 


9) 


This subroutine is not valid for METH(JD = 1,2, 5, 7. 

The finite difference method may be unstable if the bulk 
temperature method (IFM = 0) is used. The other methods 
(IFM = 1,2,3) are always stable. 

Data cards are input in groups for each control volume. 

A fluid element card (Card 1) must be input as the first 
card for each control volume. These are followed by as 
many wall element cards (Card 2) as required. The last 
wall element card for a control volume must contain a 1 in 
column 18 to indicate the end of that control volume. 

The Cp in the heat storage term is valid for a liquid flow 
or a constant pressure gas flow. For a variable pressure 


gas flow the C^ should be changed to C^. 


This can be 


approximated by the input of as C^/C^. However, in many 
cases the heat storage term will be negligible compared to 
the other terms, and can be neglected, i.e., input as 
zero. 

The expressions for fluid methods 2 and 3 are based on an 
analytical solution with negligible heat storage term. 

Fluid method 3 is for any number of wall nodes, fluid method 
2 is for only one wall node per control volume. Fluid method 
1 is valid when most of the fluid in the control is near the 
fluid outlet temperature, i.e., when hA > m C^. Fluid method 
0 is valid when the fluid in the control volume is near an 
average between the fluid inlet and outlet, i.e., when m C 
> hA. Fluid method 0 will be unstable if hA > m C 

P 

method 1, 2, and 3 are always stable. 


P 

Fluid 
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10) If no value is specified for the number of cards to follow, 
NCF0L, the subroutine assumes that 100 cards of card type 1 
will be input. This value is used to allocate core storage. 

If multiple calls to this subroutine are made within a case, 
or if the number of data cards is less than the assumed 
number of 100, the user can reduce core storage requirements 
by specifying a value for NCF0L. For the purpose of this 
estimate, three cards of card type 2 are equivalent in data 
storage requirements to one card of card type 1. For example, 
if four cards of type 1 are input and twelve cards of type 2 
are input, the value. of NCFOL is equal to (4 + 12/3) and is 
input as 8. 
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General Heat Transfer Program 



Bl 


XI 


X2 


X3 

RC0DT 

1 

■ 

FLUIDS 

r 


"S 1728 <^F ' Tp/ (A t/3600) » 


■ number of cards to follow 
fluid method 

*0 bulk temperature method 
«1 outlet temperature method 
*2 modified h method 
«3 modified ra and h method 
laminar critical Re 
turbulent critical Re 
index for laminar Nu constants 
index for turbulent Nu constants 
last card in control volume 
=0 no 
"1 yes 

inlet fluid node 
outlet fluid node 
- volume geometry option 
■-1 user specified V » XI 
= 0 rectangular V - Xl*X2*X3 

* 1 cylindrical V « -ir/4*X2^»Xl 

volume term multiplier 


element dimensions for volume 


FLUIDS - Fluid flow 


XL0 - flow development length 
B3 - flow term multiplier 

PIN - optional control volume inlet pressure 

NW - wall node 

IGA ' area geometry option 

» -1 user specified A»X4 
* 0 rectangular A*X4*X5 

® 1 cylindrical A®it*X4*X5 

B2 - area term multiplier 

- element dimensions for area 


Required curves; 
if lF/0 

f - (IF, 27,0) - friction factor 
if IFP>0 

C - (IFF, 1,0) - fluid density 
If IFP<0 

R - (|iFp 1,1,1) - gas constant 

C - (|lFP|,2,0) - fluid specific heat 
(IMD,19,1)- fluid mass flow rate 
if IH=»0 

k - (|lFP|,3,0) - fluid thermal conductivity 
u - (|lFp|,8,0) - fluid viscosity 
If IHl<0 

h - (IH,15,2) - wall heat transfer coefficient 
Note; Not valid for method - 1,2, 5, 7 



NPRNT NFMT 

22 26 


iRsningsEai 

iKMIillQHI 


74| 78 


LC 711 7I(J 
18 22 26 


IMn IC\' (Bl] [XI] [X2] 
34 38 42 46 50 


[B3] 

62 66 70 74 



I.c Nv; IH IGA [B2] [XA] [X5] 

18 22 26 30 34 38 42 46 50 54 


58 62 66 70 74 



?!l^OTl}CIBILlTy OF 'iiii 
OiaoML PAGE IS P@01 
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A. 7.0 Model Subroutine 


Section Page 

A. 7.1 0NEDMR A-91 

(one dimensional models) 
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A. 7.1 0NEDMR 


Purpose: This subroutine will compute the heat storage terms and the heat 

conduction terms for a one-dimensional thermal model composed of 
rectangular elements. 


Equation: The heat storage terms added to the heat balance are: 

0 - ifi ■> S ''i » «i' - 


where i = N1 to N2 


For each rectangular element: 

= X Y Z / 2/(N2 - Nl) for i = N1 or N2 
= X Y Z / (N2 - Nl) for i Nl or N2 
= X Y Z for i when Nl = N2 

The heat conduction terms added to the heat balance, if Nl ^ N2, 
are; . 

Q = K B (T..^^ - 


where i = Nl to (N2 - 1) . 

For a rectangular model; 

K = k Y Z / 12 AX. 

1 

where AX^ = X / (N2 - Nl) 


Data Cards: Symbolic Name 


Value 


ID 

directional option for 0NEDMC 

Nl 

node number of first element 

N2 

node number of last element 

NM 

material properties curve relative 
numbers 

X, Y, Z 

element dimensions (in) 

B 

arbitrary constant 


3 

Required Curves: p - (NM, 1, 0) - density (Ibm/ft ) 

C - (NM, 2, 0) - specific heat at constant pressure 
^ (Btu/lbm - °R) 

k - (NM, 3, 0) - thermal conductivity (Btu/hr-f t-°R) if 
Nl N2 
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1) This subroutine is available in general format 
numbers 1 and 2. 

2) Curve argument number 2 is available as: 

a) the temperature of each node for the evaluation 
of the density and specific heat 

b) the average temperature of each connected node 
pair for the evaluation of the thermal 
conductivity. 

3) The heat conduction terms are stored as two way 
terms in the heat balance. 

4) This subroutine is not valid for METHOD = 7 if 
N1 ^ N2. 

5) Care must be taken if this subroutine is used with 
the nonsequential node number option. 

Geometry Type 

Rectangular 


Subroutine 

dNEDMR 
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General Heal Tianslei Piogtam 


0NEDMR - One Dimensional Model - Rect. 


N1 - 
N2 - 
X - 







Yi 
2 
B - 


First Node 
Last Node 

Material thickness 
between nodes '-.x and 

Element Dimensions 
Constant 


Required Curves : 
p- (NM,1,0) - Density 
Cp- (NM,2,0) - Specific Heat 
N2 k - (NM,3,b) - Thermal Conductivity 

Note: May cause node number/equation 
number conflict if using nonsequential 
node number option. 



10 


18 

NPRNT 

22 

NFMT 

26 

■ 

34 

38 

42 

46 

50 

54 

58 

62 

66 

70 

74 

78 

80 

0NEDHR 


SB 



1 
















RC0DT 

1 

■ 

NCARO 
' 14 

13 

Nj. 

22 

N2 

26 

■ 

X 

3! 40 

Y 

41 50 

2 

51 60 

B 

61 70 

71 80 

1 

0NEDMR 













— 





























B 













■ 





1 













- 
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A. 8. 1 

A. 8. 2 
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A. 8.0 Miscellaneous Subroutines 


Page 

CURVES .,^95 

(curve argument) 

MOPA A-98 

(miscellaneous operations, subroutine) 
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Purpose: 


Data Cards: 


Required Curves: 


Restrictions and 


A. 8.1 CURVES 

This subroutine allows the user to store curve argument values 
from the result of a table lookup from a curve. 

Symbolic Name Value 

NREL Curve relative number, type and class 

NTYPE to be searched 

NCLAS 

NARG Curve argument number to store result 

As specified by NREL, NTYPE and NCLAS 
Notes : 

1) The subroutine CURVE is designed for general use. 

2) This subroutine must be called before any other subroutine 
which expects to use its results. 

3) The value of NARG is ignored if the type and class of the 
curve matches the atmosphere/ trajectory identifications 
shown on the data sheets. 

4) If a velocity curve is specified ( , 12, 2), the sub- 

routines will also compute and store the Mach number. 

If a Mach number curve is specified ( , 12, 3), the sub- 

routines will also compute and store the velocity. 
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FASTEMP General Heat Tianslei Piogiam CURVES - Curve Argument Loading 


- Altitude (arg. 3) NARG - Curve argument loca- Required curves: 

-,11,1 - Freestream temp. (arg. 8) tion to store value ob-^ained - (NREL.NTYPE.NCLAS) as specified 

-,11,2 - Freestream pres. (arg. 7) from curve. 

“»T2,2 - Velocity (arg.U) Note: Curve order may be important 

-,12,3 - Mach number (arg. 5) If velocity specified, Mach The curves shown will be automatically 

-,I2,t - Angle of attack (arg. 6) no. will be computed Mq^= stored in their proper curve argument. 

-,11,3 - Freestream den. (arg. 9) Voo/l*9.1\/^ 

If Mach no. specified velo- 
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Purpose : 
Equation: 


O 


A. 8. 2 M0PA 

This subroutine allows the user to perform miscellaneous 
operations for the computation and modification of dimension 
table values, curve argument values, and curve values. 

The result of the miscellaneous operation, R, is computed from 
one or two variables (VI and V2) 

R = f (VI, V2) 

depending on the operation code, I0PR. The options are: 

I0PR Result 


0 

R 

= VI 


+1 

R 

= VI + V2 


+2 

R 

= VI - V2 


+3 

R 

= VI X V2 


+4 

R 

= VI ^ V2 


+5 

R 

= minimum (Vl, 

V2) 

+6 

R 

= maximum (Vl, 

V2) 

+7 

R 

= VlV2 


+8 

R 

= average (Vl, 

V2) 



(-1.0 if Vl 

< V2 

+9 

R 

= <J 0.0 if Vl 

= V2 



(+1.0 if Vl 

> V2 

-1 

R 

= |vi 


-2 

R 

=vn 


-3 

R 

= log vl 


-4 

R 

= £n Vl 


-5 

R 

= lovi 


-6 

R 

= eVl 



Note that when the operation code is positive two variables 
are used to compute the result, otherwise only one variable is 
used. The two variables can be either dimension table values, 
curve argument values, or the result of a curve lookup. The 
result can be either a dimension table value, a curve 
argument value, or can be stored as a one dimensional curve 
value for later lookup. 

An option code, I0PT, is also provided to specify when the 
operation is performed. When the option code is zero, the 
operation is performed on all time steps. When the option 
code is +1, the operation is performed only during the first 
time step. 


A-98 

MCOOMIVELt. OOLIGI.AS COOT»>Of« ATIOM 


c-3 



FASTEMP 


Data Cards: 

Symbolic Name Value 


IRT 


IRN 


IRC 


IVIT 


IVIN 


IVIC 


I0PR 

IV2T 

IV2N 

IV2C 

I0PT 

Required Curves: 


result type 

= +1 dimension table value 
= 0 curve argunr'' .t value 

= -1 one dinensxonal curve value 
result number 

if IRT = +1, dimension table index 
if IRT = 0, curve argument number 
if IRT = -1, curve relative number 
additional result number for curve 
(IRT = -1 only) 
four digit code, ttcc, giving 
the curve type (tt) and class (cc) 
variable one type 

= +1 dimension table value 
= 0 curve argument value 
= -1 curve lookup value 
variable one number 

if IRT = +1, dimension table index 
if IRT =0, curve argument number 
if IRT = -1, curve relative number 
additional result number for curve 
(IVIT = -1 only) 
four digit code, ttcc, giving 
the curve type (tt) and class (cc) 
operation code 
variable two type (see IVIT) 
variable two number (see IVIN) 
additional variable two number (see IVIC) 
option code 


As required by the result and variable type codes. 
Restrictions and Notes: 


1) This subroutine is available only in general format number 2. 

2) This subroutine must be called before any subroutine which 
expects to use its results. 

3) This subroutines dynamically uses the dimension table for 
all variables and results, i.e., the required dimension 
table values are obtained from or stored in the dimension 
table every time step. 

4) Although this subroutine may store dimension table values, 
these values may only be referenced by this subroutine or 
other subroutines which dynamically use the dimension table. 
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R >= f 

(VI, V2) 

Type 

Value 

Operation 

Value 





+1 

dimension table value 

0 

R = VI 





0 

curve argument value 

1 

R » VI + V2 



IRT - 

result type 

-1 

curve value 

2 

R = VI - V2 



IRN - 

result number 



3 

R = VI X V2 



IRC - 

additional result 

Number 

Value 

k 

R = VI ^ V2 




number (IRT » -1 only) 

if type 

»+l dimension table index 

5 

R = min (VI, V2) 



IVIT - 

variable one type 

if type 

= 0 curve argtiment number 

6 

R = max (VI, V2) 



IVIN - 

variable one number 

if type 

““1 curve relative number 

7 

R = VlV2 



IVIC - 

additional variable one 


(also additional code gives. 

8 

R ■= av (VI, V2) 




number (IVIT = -1 only) 


U digit code, ttcc, for curve 


(-1 if VI < V2 



1V2T - 

variable two type 


type and class) 

9 

R = < 0 if VI = V2 



1V2N - 

variable two number 




(+1 if VI > V2 



IV2C - 

additional variable two 

Option 

Value 

-1 

R = |V1| 




number (1V2T ® -1 only) 

0 

performed for all time steps 

-2 

R =\Aa 



lOPR - 

operation code 

+1 

performed only during first 

-3 

R - log VI 



I0PT - 

option code 


time step 

-4 

R » In VI 







-5 

R - lOVl 







-6 

R - eVl 



INPRNT NFMT 


M0PA 

1 

0 



2 



— 








■■ 




■■■■■■ 

IHI 

■■■ 



■■■■■I 




IRT IRN IRC 
18 22 26 


IVIN IVIC I0PR IV2T IV2N IV2C I0PT 
34 38 42 46 50 54 58 


66 70 74 
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SUBROUTINE INDEX 


SUBROUTINE 

SECTION 

SUBROUTINE 

SECTION 

CURVES 

A. 8.1 

RCA 

A. 4.1 



RCC 

A. 4. 3 

DCA 

A. 2.1 

RKA 

A. 4.1 

DCB 

A. 2.1 

RKC 

A. 4. 3 

UCC 

A. 2.1 

RRA 

A. 4.1 

DCG 

A. 2. 3 

RRE 

A. 4. 2 

DCll 

A. 2. 3 



DC I 

A. 2. 3 

RRC 

A. 4. 3 

DKA 

A. 2.1 

RRD 

A. 4. 4 

UKB 

A. 2.1 

RSA 

A. 4.1 

DKC 

A. 2.1 

RSC 

A. 4. 3 

DKE 

A. 2.1 



DRA 

A. 2.1 

SCA 

A. 1.1 

DRB 

A. 2. 2 

SKA 

A. 1.1 

DRC 

A. 2. 3 

SKB 

A. 1.1 

DRD 

A. 2. 4 

SKC 

A. 1.1 

DSA 

A. 2.1 

SRA 

A. 1.1 

DSB 

A. 2.1 

SSA 

A. 1.1 

DSC 

A. 2.1 

VGA 

A. 3.1 

FLUIUB 

A. 6. 1 

VKA 

A. 3.1 



VRA 

A. 3.1 

MOPa 

A. 8. 2 

VSA 

A. 3.1 

ONHDMR 

A. 7.1 



QCCA 

A. 5.1 



QCKa 

A. 5.1 



qCRA 

A. 5.1 



QCSA 

A. 5.1 
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APPENDIX P- 

Linearization Methods 

The solution of a general set of nonlinear simultaneous equations is a diffi- 
:ult and time consuming problem. In a thermal analyzer program the heat balance 
equations must be setup and solved every time step. This is not practical even 
for a model consisting of only a few nodes. All nonlinear terms in the heat balance 
are, therefore, replaced by approximate linearized terms. These approximations are 
given below. 

(l) Radiation Approximation (all radiation subroutines) 


The radiation heat flux error may be written as 

" ^approx . ^actual 


error = 


^actual 


where 


■^actual “ 

For the forward finite difference me'-hod the temperature used is that at the 
start of a time step 

T = 

resulting in 

^actual ^ 


which is a known value and requires no approximation. For the mid finite 
difference method the temperature used is that at average time 

T = i (T° + T') 

resulting in 

'^actual 


cc (T° + T' 


If the expression is exnanded in a Taylor series about the point T° , evaluated 
at the point T', and only the linear terms in T° retained, the approximate 
radiation heat flux is given by 

q oc 3T°3 T' - 2T°^ 

approx . 

The radiation heat flux error may then be written as 


error =3(^-1) 

For the backward finite difference method the temperature used is that at 
the end of time step 
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resulting in 


q 


actual 


oc T’ 


1 | 


If this expression is expanded in a Taylor series about the point T° , 
evaluated at the point T', and only the linear terms in T° retained, the 
approximate radiation heat flux is given by 


q oci+T°3 T' - 3T°^ 

approx . 

The radiation heat fl\ix error may then be written as 

error = 4 ( — - l ) 

1^0 

Figure (0-1) shows the radiation heat flux error. It is seen that the error 
is small only when the temperature change over a time step is small. 
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T* /T^ 

FIGURE B-1 - RADIATION APPROXIMATION ERROR 
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APPENDIX C 
DATA INPUT FORMS 

Form Description Page 

Program Specification C-2 

Case Data Input C-3 

Dimension Table Input 1 C-4 

Dimension Table Input 2 C-5 

Dimension Table Input 3 ' C-6 

General Format Number 1 . C-7 

General Format Number 2 C-8 

Temporary Curve Input 1 C-9 

Temporary Curve Input 2 C-10 

Same Curve Input C-11 

Case Control Cards C-12 
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GENERAL HEAT TRANSFER 


CASE DATA INPUT 


NCASE - cose name 
NCA - number of nodes 
MPRNT - mosler print control 
N0RBIT — number of orbits 
METHf^D— finite difference method 


NSN - nonsequentiol node option 

0 ~ sequentiol nodes 

1 - nonsequentiol nodes 


3 — backward (full matrix) 

4 - backward (tridiagonal matrix) 

6 - backward (K diogonal matrix) 

8 - backward (Gauss*Seidel iteration) 


BTIME — beginning time 
FTlME - final time 

TMAX ~ moximum temperature (obsolute) 

TMIN - minimum temperoture (obsolute) 
AD02T - volue added to initiol temperotures 
TEQAL “ initiol t emperature option 

TEOAL > 0 set all to TEQAL 
TEQAL ” 0 reod in [CASE! cords) 
TEOAL ^ 0 use lost cose onswers 


I 617 


Title 


80 


CASEA ' 


I NCASE, J NCA |MPRNT^C(RBIT1METHED| IMP, I NSN I IT0 I ~ 

6 _ 8 V5 20 25| 30 35 40 45 50 55 60 _ 


CASEB 


CASEC ' 
CASET 


BTIME 

n 20 


2 , ^^'"^30 


TMAX „ 
3) 40 


41 50 


51 
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GENERAL HEAT TRANSFER 

DIMENSION TABLE INPUT 1 

HHH 

im 


DIMENS 




1 6 

KE91 

kbEI 

mmmmm 

26 40 

41 55 

56 70 

VALUES 

01 







05 







09 





■ 


13 





s 


17 







21 








25 







29 





■ 


33 





B 


37 





1 


41 





— 1 


45 







49 







53 







57 







61 





~i 


65 







69 







73 







77 







81 





■ 

■ 

85 
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GENERAL HEAT TRANSFER 

DIMENSION TABLE INPUT 3 

1 6 

NO. 

8-10 

11 25 

26 40 

41 55 

56 70 

VALUES 

189 





B 

m 

193 







197 







201 







205 







209 







213 







217 







221 







225 







229 







233 







237 







241 







245 







249 







253 







257 







261 







265 







269 







273 







277 







281 







285 






C-6 

tkfCDOIVIVeLt. OOUGL.A9 CO ft r* on A TIOmi 




































FASTEMP 


REPORT MOC A3791 



c-" 

MCootviveuL ooc/GLAa coupon A noun 











FASTEMP 


REPORT MDC A3791 


FASTEMP 

GENERAL HEAT TRANSFER PROGRAM 

GENERAL FORMAT NUMBER 2 


10 


Bai 

m 

m 


m 

m 

m 



H 

m 

n 


m 

m 

m 











































































































■ 




















B 





















' 






















































1 


















— 



— 




















- 






m 




















B 








































1 

















































































m 




















B 





























m 




















B 
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GENERAL HEAT TRANSFER 

TEMPORARY CURVE INPUT 1 

1 6 


mjg 


22 






TCIDNT 










16 7 80 

TCTITL 


11 25 26 40 41 55 56 70 














■ 



1 6 

la 


8M 

fm 


NDIMN 

30 

NPT|1 

NARG 1 
38 


TCIDNT 


Hi 









1 6 7 


TCTITL 


11 25 26 40 41 55 56 70 

















■■ 




HR 


NDIMN 

30 




•CIDNT 










1 6 7 80 

TCTITL 


II 25 26 40 41 55 56 70 
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1 

GENERAL HEAT TRANSFER 

TEMPORARY CURVE INPUT 2 



NTTYPE 

14 

14CLAS 

18 

22 


i 

■EQ 



H 


n 


VFTS4 

58 


m 

sn 

!■■■■ 

TCIDNT 










r 








16 7 80 

TCTITL 


11 25 26 40 41 55 56 70 























































1 




.. 
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